HD-THEP-02-04 
hep-ph/0202012 



S'-Matrix Unitarity, Impact Parameter Profiles. 
Gluon Saturation and High-Energy Scattering 

A. I. Shoshi 1 ^, F. D. Steffen 1 ^, and H. J. Pirner 1 ' 2 ^ 

1 Institut fur Theoretische Physik, Universitdt Heidelberg, 
Philosophenweg 16 & 19, D-69120 Heidelberg, Germany 

2 Max-Planck-Institut fur Kernphysik, Postfach 103980, 
D-69029 Heidelberg, Germany 

Abstract 

A model combining perturbative and non-perturbative QCD is developed 
to compute high-energy reactions of hadrons and photons and to investigate 
saturation effects that manifest the S-matrix unitarity. Following a functional 
integral approach, the S'-matrix factorizes into light-cone wave functions and 
the universal amplitude for the scattering of two color-dipoles which are rep- 
resented by Wegner- Wilson loops. In the framework of the non-perturbative 
stochastic vacuum model of QCD supplemented by perturbative gluon ex- 
change, the loop-loop correlation is calculated and related to lattice QCD 
investigations. With a universal energy dependence motivated by the two- 
pomeron (soft + hard) picture that respects the unitarity condition in impact 
parameter space, a unified description of pp, np, Kp, j*p, and 77 reactions is 
achieved in good agreement with experimental data for cross sections, slope 
parameters, and structure functions. Impact parameter profiles for pp and 
^* L p reactions and the gluon distribution of the proton xG(x, Q 2 , \ b±\) are cal- 
culated and found to saturate in accordance with S'-matrix unitarity. The 
cm. energies and Bjorken x at which saturation sets in are determined. 
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1 Introduction 



One of the challenges in quantum chromodynamics (QCD) is the description and 
understanding of hadronic high-energy scattering. Since the momentum transfers 
can be small, the QCD coupling constant is too large for a reliable perturbative 
treatment. Non-perturbative QCD is needed which is rigorously only available as 
a computer simulation on Euclidean lattices. Since lattice QCD is limited to Eu- 
clidean space-time, it cannot be applied in Minkowski space-time to describe par- 
ticles moving near the light-cone. Only static properties of hadrons such as the 
hadron spectrum or the phenomenon of confinement can be accessed within lattice 
QCD until now. 

An interesting phenomenon observed in hadronic high-energy scattering is the 
rise of the total cross sections with increasing cm. energy. While the rise is slow in 
hadronic reactions of large particles such as protons, pions, kaons, or real photons 
it is steep if only one small particle is involved such as an incoming virtual photon p], 
|3j or an outgoing charmonium ||. 

This energy behavior is best displayed in the proton structure function F 2 (x, Q 2 ) 
that is equivalent to the total 7*p cross section. With increasing photon virtuality 
Q 2 , the increase of F 2 (x, Q 2 ) towards small Bjorken x becomes significantly stronger. 
Together with the steep rise of the gluon distribution in the proton xG(x, Q 2 ) with 
decreasing x, the rise of the structure function F 2 (x,Q 2 ) towards small x |2|,|| is 
one of the most exciting results of HERA. 

Postulating the unitarity of the S- matrix, SS^ — S'S = H, asymptotic lim- 
its on the growth of total hadronic cross sections have been derived such as the 
Froissart-Lukaszuk-Martin bound ||. This limit allows at most a logarithmic en- 
ergy dependence at asymptotic energies. Analogously, the rise of the total 7*p cross 
section is expected to slow down due to parton saturation effects reflecting S'-matrix 
unitarity. In fact, it is a key issue to determine the energy at which unitarity limits 
become significant. 

A phenomenologically very successful and economical description of the energy 
dependence in both hadron-hadron reactions and 7*p reactions is given by the two- 
pomeron model of Donnachie and Landshoff ||. In this picture, the energy depen- 
dence of the cross sections at high energies results from the exchange of a soft and a 
hard pomeron, the first of which dominates in hadron-hadron and 7*p reactions at 
low Q 2 and the second one in 7*p reactions at high Q 2 || . Both pomerons carry 
by definition the quantum numbers of the vacuum and may be related to a glueball 
trajectory || or a gluon ladder ||. The two-pomeron model, however, explicitly vi- 
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olates the Froissart-Lukaszuk-Martin bound || at asymptotic energies and does not 
contain parton saturation. A model motivated by the concept of parton saturation 
is the one of Golec-Biernat and Wiisthoff || which allows very successful fits to 7*p 
data but cannot be applied to hadron-hadron reactions. 

In this work, we develop a model combining perturbative and non-perturbative 
QCD to compute high-energy reactions of hadrons and photons with special empha- 
sis on saturation effects that manifest the S-matrix unitarity. Aiming at a unified 
description of hadron-hadron, photon-hadron, and photon-photon reactions involv- 
ing real and virtual photons as well, we follow the functional integral approach to 
high-energy scattering in the eikonal approximation |10|-[T3|, in which the S'-matrix 
element factorizes into the universal correlation of two light-like Wegner- Wilson 
loops Sdd and reaction-specific light-cone wave functions. The light-like Wegner- 
Wilson loops describe color-dipoles given by the quark and antiquark in the meson 
or photon and in a simplified picture by a quark and diquark in the baryon. Conse- 
quently, hadrons and photons are described as color-dipoles with size and orientation 
determined by appropriate light-cone wave functions [T^,[T3|. Thus, the loop-loop 
correlation function Sdd is the basis for our unified description. 

We evaluate the loop-loop correlation function Sdd in the approach of Berger 
and Nachtmann |l4j . In this approach, the S'-matrix unitarity condition is respected 



as a consequence of a matrix cumulant expansion and the Gaussian approximation 
of the functional integrals. We explicitly review the Berger-Nachtmann approach as 
it is crucial for our loop-loop correlation model and our investigation of saturation 
effects. 

We express the loop-loop correlation function Sdd in terms of the gauge-invariant 
bilocal gluon field strength correlator integrated over two connected surfaces. These 
surfaces enter from an application of the non-Abelian Stokes' theorem, in which 
the line integrals are transformed into integrals over surfaces bounded by the loops. 
We use for the first time explicitly minimal surfaces. This surface choice is usually 
used to obtain Wilson's area law in Euclidean space I5|, 16|. The simplicity of 



the minimal surfaces is appealing. It allows us to present the explicit computation 
of Sdd compactly in this work and to extract an analytic structure of the non- 



perturbative contribution to the dipole-dipole cross section in [17 



We decompose the gluon field strength correlator into a non-perturbative and a 
perturbative component. The stochastic vacuum model (SVM) |15 is used for the 



non-perturbative low frequency background field and perturbative gluon exchange for 
the additional high frequency contributions. This combination allows us to describe 
long and short distance correlations in agreement with Euclidean lattice calculations 



of the gluon field strength correlator |T8|,|T^]. Moreover, if applied with the minimal 



2 



surface in Euclidean space-time, this two component ansatz leads to the static quark- 
antiquark potential with color-Coulomb behavior at short distances and confining 
linear rise at long distances jnj . In this way, a connection of high-energy scattering 
to lattice simulations of QCD and the QCD string tension is established. 

We use in the non-perturbative component the exponential correlation function 
directly from lattice QCD investigations of the correlator This correlation 

function stays positive for all Euclidean distances and, thus, is compatible with a 



spectral representation of the correlation function |2(| . This means a conceptual im- 



provement since the correlation function that has been used in earlier applications 



of the S VM becomes negative at large distances |L2], [14], pTH25|| . Besides the concep- 
tual and computational advantages, the new combination — exponential correlation 
function and minimal surfaces — allows a successful phenomenological description 
of the slope parameter B(s), the differential elastic cross section da el /dt(s,t), and 
the elastic cross section a el (s) as shown in this work. 

We introduce the energy dependence into the loop-loop correlation function Sbd 
in order to describe simultaneously the energy behavior in hadron-hadron, photon- 
hadron, and photon-photon reactions involving real and virtual photons as well. 
Motivated by the two-pomeron picture of Donnachie and Landshoff ||, we ascribe 
to the soft and hard component a weak and strong energy dependence, respectively. 
Including multiple gluonic interactions, we obtain an S-matrix element with a uni- 
versal energy dependence that respects unitarity in impact parameter space. 

We adjust the model parameters to reproduce a wealth of high-energy scattering 
data, i.e. total, differential, and elastic cross sections, structure functions, and slope 
parameters for many different reactions over a large range of cm. energies. In this 
way, we have confidence in our model predictions for future experiments (LHC, 
THERA) and for energies beyond the experimentally accessible range. 

To study saturation effects that manifest the S'-matrix unitarity, we consider 
the scattering amplitudes in impact parameter space, where the S'-matrix unitar- 
ity imposes rigid limits on the impact parameter profiles such as the black disc 
limit. Having confirmed that our model respects the unitarity condition in impact 
parameter space, we compute profile functions for proton-proton and longitudinal 
photon-proton scattering that explicitly saturate at the black disc limit at high en- 
ergies. These profiles provide also an intuitive geometrical picture for the energy 
dependence of the cross sections. 

Using a leading twist, next-to-leading order DGLAP relation, we estimate the 
impact parameter dependent gluon distribution of the proton xG(x, Q 2 , \b±\) from the 
profile function for longitudinal photon-proton scattering. We find low-x saturation 
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of xG(x, Q 2 , \b±\) as a manifestation of the S- matrix unitarity. The implications on 
the integrated gluon distribution xG(x, Q 2 ) are studied and compared with comple- 
mentary investigations of gluon saturation. 

With the profile function saturation and the intuitive geometrical picture gained 
in impact parameter space, we turn to experimental observables to localize satu- 
ration effects in the cross sections and to interpret the energy dependence of the 
cross sections. We compare the model results with the experimental data and pro- 
vide predictions for future experiments and saturation effects. Total cross sections 
a tot , the structure function of the proton F2, slope parameters B, differential elastic 
cross sections do~ el /dt, elastic cross sections a el , and the ratios a el /a tot and a tot / B are 
considered for proton-proton, pion-proton, kaon-proton, photon-proton, and photon- 
photon reactions involving real and virtual photons as well. 

The outline of the paper is as follows: In Sec. 0, the model is developed and 
the model parameters are given. Going to impact parameter space in Sec. [|, we 
confirm the unitarity condition in our model and study the impact parameter profiles 
for proton-proton and photon-proton scattering. The impact parameter dependent 
gluon distribution of the proton xG(x,Q 2 , \b±\) and gluon saturation are discussed 
in Sec. |j. Finally, in Sec. |5], we present the phenomenological performance of the 
model and the saturation effects in the experimental observables. The appendices 
present explicitly the used hadron and photon light-cone wave functions and the 
analytic continuation of the non-perturbative correlation functions from Euclidean 
to Minkowski space-time. 



2 The Loop-Loop Correlation Model 



The T-matrix is the central quantity in scattering processes. It enters every ob- 
servable we intend to look at and is obtained from the ^-matrix by subtracting the 
trivial case in which the final state equals the initial state, 



% = 5 fi + i(2kY5\P s - P t )T f 



(2.1) 



where Pi and Pf represent the sum of incoming and outgoing momenta, respectively. 
We compute the T-matrix in a functional integral approach developed for parton- 



parton scattering []10[] in the eikonal approximation and extended to gauge-invariant 



loop-loop scattering |11-13|. In this approach, the T-matrix element for the reaction 
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ab — > cd factorizes as follows 



T ab ^ cd (s,t) = 2is Jd 2 b ± e^ Jdz l d 2 r 1 Jdz 2 d 2 r 2 

x ip*(zi,rx) ip* d (z 2 ,r 2 ) 1 - S DD (b±,z 1 ,f 1 ,z 2 ,f 2 ) ipa^,^) ip b (z 2 ,r 2 ) , (2.2) 



where the loop-loop correlation function 

S DD (b±,z 1 ,r 1 ,z 2 ,r 2 ) 



(w[c x ]w[c 2 ) 



G 



(2.3) 



describes the elastic scattering of two color-dipoles (DD) with transverse size and 
orientation r*j and longitudinal quark momentum fraction Z{ at impact parameter 
b±, transverse momentum transfer q± (t = —q 2 ) and cm. energy squared s. In this 
framework, the color-dipoles are given by the quark and antiquark in the meson or 
photon and in a simplified picture by a quark and diquark in the baryon. Conse- 
quently, the hadrons and photons are characterized by the light-cone wave functions 
ipafi and ip C) d that describe the r*j and ^ distribution of the color-dipoles. Concen- 
trating in this work on reactions with a = c and b = d, only squared wave func- 
tions |V>i(^i,ri)| 2 := if)*(zx,n) i> a (zi,ri) and \^ 2 (z 2 ,r 2 )\ 2 := ^ d {z 2 ,r 2 )^ h {z 2 ,r 2 ) are 
needed. We use for hadrons the phenomenological Gaussian wave function [^5|,|26|| 
and for photons the perturbatively derived wave function with running quark masses 
mf(Q 2 ) to account for the non-perturbative region of low photon virtuality Q 2 
as discussed explicitly in Appendix |A|. 

The path of each color-dipole is represented by a light-like QCD Wegner- Wilson 
loop [^J 

1 



]22 



W[C h2 ] = — TrPexp 



-ig (p dz^G^z) 

Cl,2 



(2.4) 



where N c is the number of colors, Tr the trace in color space, g the strong cou- 
pling, and Gix{z) = Q^{z)t a the gluon field with the SU(N C ) group generators t a 
that demand the path ordering indicated by V. Quark-antiquark dipolesQ are rep- 
resented by loops in the fundamental SU(N C = 3) representation. In the eikonal 
approximation to high-energy scattering the q and q paths form straight light-like 
trajectories. Figure illustrates the space-time (a) and transversal (b) arrangement 
of these loops. The world line C\ (C 2 ) is characterized by its light-cone coordinate 
x~ = x° — x 3 = (x + = x° + x 3 = 0), the transverse size and orientation f\ (f 2 ) 
and the longitudinal quark momentum fraction Zi(z 2 ) of the corresponding dipole. 
The impact parameter between the loops is 



n q + (1 - z x )rx - f 2q - (1 - z 2 )r 2 = r lcm - r 2cm 



(2.5) 



*or equivalcntly quark-diquark (q — qq) or antiquark-diantiquark systems (q — qq) 
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as shown in Fig. |I]b, where f q (fig) is the transverse position of the quark (antiquark) 
in loop i, fi = f iq - f ig , and f icm = Zif iq + (1 - z^fg. 

The QCD vacuum expectation value (. . .)g in the loop-loop correlation function 
( [2.3|) represents functional integrals [|13[ in which the functional integration over the 



fermion fields has already been carried out as indicated by the subscript G. The 
model we use for the QCD vacuum (see Sec. |2.2|) describes only gluon dynamics 
and, thus, implies the quenched approximation that does not allow string breaking 
through dynamical quark-antiquark production.^ 



2.1 The Loop-Loop Correlation Function 



To compute the loop-loop correlation function ( |2.3|) , we transform the line integrals 
over the loops into integrals over surfaces S\ >2 with dSip = by applying 



the non-Abelian Stokes' theorem [fT3|, |28 

'w[Ci]w[c 2 



= (^rTrP s exp 

G \N r 



x — TrP 5 exp 



i| ! da^{x x )Ql v {o x ,x x -C Xl0x )t a 
i| / do- pa (x 2 )G b pa (o2,x 2 ;C X202 )t 

1 J So. 



(2.6) 



where the gluon field strength tensors, Q^ v {x) = G^ v (x)t a , are parallel transported 
to the point o along the path C xo 



G^{o, x; C xo ) = $(x, o; C xo ) l Q lxv {x)^{x, o; C xo ) 
with the QCD Schwinger string 



$(x,o;C xo ) = Pexp 



Cxo 



(2.7) 



(2- 



In ( p.6|) , Vs indicates surface ordering and 0\ and o 2 are the reference points on 
the surfaces Si and S 2 , respectively, that enter through the non-Abelian Stokes' 
theorem. In order to ensure gauge invariance in our model, the gluon field strengths 
associated with the loops must be compared at one reference point o. Therefore, we 
require the surfaces S\ and 5*2 to touch at a common reference point 01 = 02 = 0. 



Following the approach of Berger and Nachtmann |TJj], the product of the two 
traces (Tr) over iV c x iV c matrices in (|2.6|) is expressed as one trace (Tr 2 ) that acts 



5 The quenched approximation becomes explicit in the linear rise of the dipole-proton and dipole- 
dipole cross-section with growing dipole size obtained in our model. 
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Figure 1: Space-time (a) and transverse (b) arrangement of the Wegner- Wilson loops. 
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in the ^-dimensional tensor product space of two fundamental SU(N C ) represen- 
tations 



{w[C x \W[C^ q = (^Tr 2 



P 5 exp[-*! / ,la""U i )g' l ' ll Au..r L :C. l . i „)t"} : 



X 



.9 



Using the identities 



11 <g> P 5 exp[-z| / ^ CT (x 2 )^ (T (o,x 2 ;C a;2 o)t 6 



exp ( t a ) <g> 11 = exp ( t a ® 11 ) 
11 (g) exp ( f 1 ) = exp ( 11 <g> t a ) 



G 



• (2-9) 



(2.10) 
(2.11) 



the tensor products can be shifted into the exponents. With the matrix multiplica- 
tion in the tensor product space 



a ib 



(t a ® 11) (t b <g> 11) = t a t 



11 



t° <g> t 6 



(2.12) 



(t a ® 11) (11 ® t b ) 
and the vanishing commutator 

[t a ® 11, 11 ® t 6 ] = , (2.13) 
the two exponentials in (|2.9|) commute and can be written as one exponential 



(W[C 1 ]W[C 2 ]) g = ^Tr 2 V s ex V 



-i| I da^(x)G^(o,x;C x 



G 



(2.14) 



with the following gluon field strength tensor acting in the iVf- dimensional tensor 
product space 



(o, x; C x 

O 



g%(o,x-C xo )(t a ® 11) for i G S! 
(o^C^Xll <g> t fe ) for x G S 2 



(2.15) 



In ( p. 14 ), the surface integrals over Si and S* 2 are written as one integral over the 
combined surface S = Si + SV For the evaluation of ( |2.14| ), the linearity of the 
functional integral, (Tr ■ - ■) = Tr(- ■ •), and a matrix cumulant expansion is used as 
explained in (cf. also P|) 



P 5 exp -z- J da(x)G(o,x;C xo ) 

00 i r 

= exp[ ^?( _i 2^ n / da ( Xl ) ' ' ' K r 



(x^ , ■ ■ ■ , x n )] 



(2.16) 



n=l 



where Lorentz indices are suppressed to lighten notation. The cumulants K n consist 
of expectation values of ordered products of the non- commuting matrices Q(o, z\ C zo ). 
The leading matrix cumulants are 

K^x) = (g(o,x;C x )) G , (2.17) 
K 2 {x 1 ,x 2 ) = (VsG(o,x 1 ;C xl )Q(o,x 2 ;C X2 )) G 

-^(g(o,x 1 ;C Xl )) G (g(o,x 2 ;C X2 )) G + (1^2)^) . (2.18) 

Since the vacuum does not prefer a specific color direction, K\ vanishes and K 2 
becomes 

K 2 (x l7 x 2 ) = (V s G(o, xi, C Xl )G(o, x 2 ; C X2 )) G . (2.19) 

Now, we restrict the functional integral associated with the expectation values (. . .) G 
to be a Gaussian functional integral. Consequently, all higher cumulants, K n with 
n > 2, vanishF] and the loop- loop correlation function can be expressed in terms of 
K 2 



W[C X \W[Q 



G 



^Tr 2 exp 



^- \ do^ v [xi) I da pa (x 2 )(Vs G^u{o, xi] C Xl0 )G p<7 (o, x 2 , C X2D 
° Js Js x 1 G 



[2.20) 



Using definition ( 2.15|) and the relations (|2.12|) , we now redivide the exponent in 
( p.20| ) into integrals of the ordinary parallel transported gluon field strengths over 
the separate surfaces Si and S 2 



W[C 1 ]W[C 2 ]) g = ^Tr 2 exp 



9 2 



da^( Xl ) da^x 2 )P s (g^o, Xl ;C xlo )g%(o,x 2 ;C^ )) (t a ® t b ) 
Js 2 x 1 G 

^ / da^( Xl ) [ da^ix^Vsig^x^C^g^x^C^)) (t a ® t b ) 
° Js 2 J St x ' 1 G 

^ / da^{x x ) [ da^(x 2 )V s (g; v (o,x 1 -,C Xl0 )g b pa (o,x 2 ,C X20 )) (t a t b ® 11) 

° J Si J Si x 1 G 



9 2 



IS! JS! 

da^(x x ) [ da^(x 2 )V s (g;Ao,x 1 ;C Xl0 )g b (T (o,x 2 ;C X20 )) (11 ® t a t b 
s 2 Js 2 x 1 G 



(2.21) 



6 We are going to use the cumulant expansion in the Gaussian approximation also for perturba- 
tive gluon exchange. Here certainly the higher cumulants are non-zero. 
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Due to the color-neutrality of the vacuum, the gauge- invariant bilocal gluon field 
strength correlator contains a 5-function in color-space, 

,2 



/ 9 



1 



\^^u(o,x 1 ;C Xl0 )g b pa (o,x 2 ;C X20 )) G =: -^(n, x 2 , o; C X10 , C X20 ) , (2.22) 

which makes the surface ordering Vs in (|2.21|) irrelevant. The quantity F^ vp(7 will 
be specified below. With ansatz (|2.22j) and the definition 



XSiS s '■= / da^(xi) / da pa (x 2 )F fllJpa (x 1 , x 2 , o; C XlD , C X2Q ) , (2.23) 

Eq. ( ggTJ ) reads 



G 



— i 



(Xs lS2 + Xs 2Sl ) (t a ® t°) 



+ Xs lSl (t a t a <g> 11) + X5 2 5 2 (U ® n a ) 



(2.24) 



Our ansatz for the tensor structure of F, 



/ivpa 



see (|273T| ), ( ggg ), and ( |2~37D 



leads to XSiSi — Xs 2 s 2 = f° r light-like loops, as explained in Sec. |2T3| , and also to 
XSiS 2 = Xs 2 S! ='■ X- F° r the evaluation of the trace of the remaining exponential, 
we employ the projectors 



(-Pa)(a 1 a 2 )(/3 1 /3 2 ) = 2 (^ift ^*2/3 2 _ ^Qi/3 2 ^a 2 /3i ) , 



(2.25) 
(2.26) 



that decompose the direct product space of two fundamental SU(N C ) representa- 
tions, in short N c , into the irreducible representations 



N c ® N c = (N c + l)N c /2 © N C (N C - l)/2 . 



With the identity 



t a ® t a 
and the projector properties 



2iV, 



2AA, 



(2.27) 
(2.28) 



P = P 

1 s.a 1 s,a 1 



Tr 2 P S = (N C + I) N c /2 , and Tr 2 P a = (N c - l)N c /2 , (2.29) 



we find for the loop-loop correlation function in the fundamental SU(N C ) represen- 
tation 



W[Ci]W[C 2 ] 



a 



N c + 1 
2N r 



exp 



.Nr. 



—I- 



2N r 



-X 



N 



2N r 



exp 



,N C + 1 
2N r 



-X 



[2.30) 



and recover, of course, for N r = 3 the result from 14 
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2.2 Perturbative and Non-Perturbative QCD Components 



We decompose the gauge-invariant bilocal gluon field strength correlator ( |2.22| ) into 
a perturbative (P) and non-perturbative (NP) component 



f np , f p 

fii/pa 1 fiupcr 



(2.31) 



^ pupa 

Here, i 7 ^,,. gives the low frequency background field contribution modelled by the 
non-perturbative stochastic vacuum model (SVM) @ and F p pa the additional high 
frequency contributions described by perturbative gluon exchange. Such a decom- 
position is supported by lattice QCD computations of the Euclidean field strength 
correlator plffSfl. 



In the SVM, one makes the approximation that the correlator Fj^ depends 
only on the difference z := X\ — x 2 but not on the reference point o and the curves 
C Xl0 and C X2Q |fT5| . Then, the most general form of the correlator that respects 
translational, Lorentz, and parity invariance reads in four-dimensional Minkowski 



space-time 11,12 



Fpvpa( Z ) 



YZTJ) ° 2 { K (dw9™ ~ 9iut9vp) D(z 2 /c 



3(iV c 2 _ 



5 



3(JV* 

-a- 



G 9 



\ z o9p.p 

d A k 



-ikz 



K (fl^SW -9na9up) D(k 2 



D~'J (2vr) 4 

kuk a 9pp kykpg^u + k^kpg ua kpk a g up D^(k / 



. (2.32) 



Here, a is the correlation length, G 2 : = (^^^(0)^^(0)) is the g/uon conden- 



sate |p0| , k determines the non-Abelian character of the correlator, D and -Di are 
correlation functions in four dimensional Minkowski space-time, and 



tfcz 



(2.33) 



In the case of k ^ 0, the Euclidean version of Fp,J}a(z) in (|2.32j ) leads to confine- 
ment and does not fulfill the Bianchi identity. In contrast, the Euclidean version of 
Fp U p { a C) (z) fulfills the Bianchi identity but does not lead to confinement ||I5| . There- 
fore, we call the tensor structure multiplied by k non-Abelian or confining (c) and 
the one multiplied by (1 — k) Abelian or non-confining (nc). 
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The non-perturbative correlator was originally constructed in Euclidean space- 



time ||15|| . The transition to Minkowski space-time is performed by the substitution 
5 — > —Qp P and the analytic continuation of the Euclidean correlation functions 
to real time, D E —> D and Df —>■ D\ [11, 12]. Euclidean correlation functions are 



accessible together with the Euclidean correlator in lattice QCD dig, [L9|[. We adopt 
for our calculations the simple exponential correlation functions specified in four 
dimensional Euclidean space-time 



D*(Z l /a z ) = D?(Z'/a z ) = exp(-\Z\/a) , (2.34) 

that are motivated by lattice QCD measurements of the gluon field strength corre- 
lator ||T%|, . These correlation functions stay positive for all Euclidean distances 
Z. In earlier applications of the SVM, a different correlation function D E has been 
used that becomes negative at large distances jl2|, |i~4l, |21~H2"5| . Such a negative part 



is not compatible with a spectral representation of the correlation function [20 



By analytic continuation of ( |2.34|) we obtain the Minkowski correlation functions 
in (|2.32|) as shown in Appendix M. 



Treating the vacuum fluctuations as a Gaussian random process, the non-perturbative 
Euclidean correlator leads to the following explicit expression for the QCD string 
tension |15] 

" = ^rl ^W) = ^~ (2.35) 

with the exponential correlation function ( |2.34|) used in the final step. The QCD 
string tension a characterizes the confining quark- ant iquark potential and can be 
computed from first principles within lattice QCD j3l|. Thus, relation ( |2.35| ) puts an 
important constraint on the three fundamental parameters of the non-perturbative 
QCD vacuum — a, and k — and eliminates one degree of freedom. 



While a non-perturbative model must be used to describe the low frequency 
contributions, the perturbative component F^ vpo . is computed from the gluon prop- 
agator in Feynman-'t Hooft gauge 

with an effective gluon mass introduced to limit the range of the perturbative 
interaction in the infrared (IR) region. 

In leading order in the strong coupling g, the bilocal gluon field strength correla- 
tor is gauge-invariant already without the parallel transport to a common reference 
point so that F^ upa depends only on the difference z := X\ —x 2 - In this order, 0(g 2 ), 
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we obtain 



£1 

7T 2 2 



5 



3 



Dp(z< 



g 2 f d A k e _ ikz 



7T 2 J (2tt)< 



kyk^g^p kykpg^u + k^k p g vtJ k^k a g V p Dp(k ) 

(2.37) 



with the perturbative correlation function 

d 
dk 2 



D' p (k 2 ) := I d 4 zD P (z 2 )e ikz 



k 2 — uiq 



(2.38) 



The tensor structure in (2.37 ) is identical to the non-confining tensor structure 
in the non-perturbative component ( |2.32| ). Together with the perturbative correla- 
tion function in Euclidean space-time, it leads to the non-confining color-Coulomb 
potential that is dominant for small quark- ant iquark separations 



In the final step of the computation of % hi the next section, the constant coupling 



g is replaced by the running coupling 



47ra s (il) 



12tt 



(33 — 2iV/) In [(|4| 



-2 



M 2 )/A 2 QCD ] 



[2.39) 



with the renormalization scale provided by \z]_ | that represents the spatial separation 
of the interacting dipoles in transverse space.Q In ( |2.39| ), Nf denotes the number of 
dynamical quark flavors, which is set to Nf = in agreement with the quenched 
approximation, Aqcd = 0.25 GeV, and M 2 allows us to freeze g 2 for \z±\ — > oo. 
Relying on a low energy theorem [^3[, we freeze g 2 at the value at which the results 
for the potential and the total flux tube energy of a static quark-antiquark pair 



coincide in our model [16] 



2.3 Evaluation of the %-Function with Minimal Surfaces 



For the computation of the x-function ( |2.23|) 



X • Xc "X-nc ' "X. 



*-nc 

r 2 



-i^-f da^{x x ) I da^(x 2 ) (fS? + FZt c) + F?} , (2.40) 

J Si J S2 



7 Timc-like or light-like separations do not appear in the final expression for \- They are 



integrated out as explained in Sec. 2.3 
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one has to specify surfaces S± t2 with the restriction 95*1,2 = Ci,2 according to the 
non-Abelian Stokes' theorem. As illustrated in Fig. [I], we put the reference point o at 
the origin of the coordinate system and choose for Si }2 the minimal surfaces that are 
built from the areas spanned by the corresponding loops Ci t2 and the infmitesimally 
thin tube which connects the two surfaces S\ and 5*2. Since the tube contributions 
cancel mutually, this choice makes the calculation explicitly independent of the 
reference point o and of the paths C XlQ and C X2Q . 

The minimal surfaces S\ and £2 shown in Fig. [l| can be parametrized with the 
upper (lower) subscripts and signs referring to Si (S 2 ) as follows 



Si 

(2) 

where 



(e) 



x£ (u, v) 

(2) 



( 1 ^ 



V A 1 / 



(2«J) 



(e) 



(2) 



/ \ 



'lq 

(2q) 







V u / 



ue[-T,T], ve [0,1] 



/ \ 



and 



n 

(2) 



C2) 





V u / 



(2.41) 



^2.42) 



The infmitesimally thin tube is neglected since it does not contribute to the \- 
function as already mentioned. The computation of the x-function requires only the 
parametrized parts of the minimal surfaces ( p. 41 ), the corresponding infinitesimal 
surface elements 



dx^ dx u dx^ dx u 



du dv dv du 



du dv 



rC ri - rr rim 



(e) (2) 



(2) (e) 



du dv 



(2.43) 



and the limit T — > 00 which is appropriate since the correlation length a is much 
smaller (see Sec. |2.5|) than the longitudinal extension of the loops. 



Starting with the confining component 

^2 



X 



NP 



7T 

T 



da^{x 



Si 



da prT (x 2 ) F, 



[11/ pa I 



Xi - X 2 ) 



S-2 



7C 2 G 2 K 

12(iV| - i; 



Si 



da^ixx) / da pa {x 2 ){g w g^-g^g vp )iD{z 2 /a 2 ) , (2.44) 



S-2 



one exploits the anti-symmetry of the surface elements, da pv = —da up , and applies 
the surface parametrization ( 2.41 ) with the corresponding surface elements ( |2.43| ) 
to obtain 



xT 



TT 2 GoK 



3(JV, 



2 (ri • f 2 ) 



1 pi pT pT 

dv\ / dv 2 lim / du\ / du 2 iD(z 2 / a 2 ) , (2.45) 

J Q T^oo_/_ T J_ T 
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where 

z^ = x± — x% = uxn 1 ^ — u 2 n,Q + r^ q — r% q + V\r^ — v 2 r% , (2.46) 

and the identities n 9 ■ r 2 = r\ ■ n Q = and n ffi ■ n e = 2, evident from ( |2.42j ), have 
been used. Next, one Fourier transforms the correlation function and performs the 
U\ and u 2 integrations in the limit T — > oo 



f T f T 
lim / du\ / du 2 iD(z 2 /a 2 ) 

T^cvJ-T J-T 

r d^k ~ r T r T 

= / 7rm^( fc2 ) lim / du i / du 2 e~ ihz 

J (2vr) 4 't^J_ t J_ t 

= J -j0- ib(k 2 ) expf-i^K, - r% + - v 2 r%)] 5(k° - k 3 ) 5(k° + k 3 ) 

= - z\D< 2 > (f lg + vin - r 2q - v 2 f 2 ) , (2.47) 

where iD^ is the confining correlation function in the two-dimensional transverse 
space (cf. Appendix |B|) 

^ (2)(4) = / ^^ &2 \k) ■ (2.48) 

The contributions along the light-cone coordinates have been integrated out so that 
Xc F is completely determined by the transverse projection of the minimal surfaces. 
Inserting (|2~47D into (|2~45p , one finally obtains 



ir 2 G f 1 f 1 



With L)( 2 \k±) obtained from the exponential correlation function (|2.34 ), cf. Ap- 
pendix ||, we find 

iD {2 \z L ) = 27ra 2 [l + (|4|/a)]exp(-|4|/a) (2.50) 

which is positive for all transverse distances. 

As evident from the v% and v 2 integrations in (|2.49 ) and Fig. |l]b, there are contri- 
butions from the transverse projections of the minimal surfaces (Si >2 )± connecting 
the quark and antiquark in each of the two dipoles. We interpret these contribu- 
tions as a manifestation of the strings that confine the quarks and antiquarks in the 
dipoles and understand, therefore, the confining component \c P as a string-string 
interaction. This component gives the main contribution to the scattering amplitude 



in the non-perturbative region [17 
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Due to the truncation of the cumulant expansion or, equivalently, the Gaussian 
approximation, a considerable dependence of \c P on the specific surface choice is 
observed. In fact, a different and more complicated result for Xc F was obtained with 
the pyramid mantle choice for the surfaces Si 2 m earlier applications of the SVM 



to high-energy scattering [|l2j, |14], p|-p5| . However, we use minimal surfaces in line 
with model applications in Euclidean space-time: If one considers the potential of a 
static quark-antiquark pair, usually the minimal surface is used to obtain Wilson's 
area law |l5|,|T6fl. Moreover, the simplicity of the minimal surfaces allows us to 
give an analytic expression for the leading term of the non-perturbative dipole- 
dipole cross section [T7|. Phenomenologically, in comparison with pyramid mantles, 
the description of the slope parameter B(s), the differential elastic cross section 
da el I dt(s , t) , and the elastic cross section a el (s) can be improved with minimal 
surfaces as shown in Sec. [|. 

Continuing with the computation of the non-confining component 
Xnc -=-^ I da^{x x ) [ da pr7 (x 2 ) F^ c \z = x x - x 2 ) 



1 Si 

Wf^i L n) lr™ 



(2.51) 



x 



d 4 k 



-ikz 



kykfjQ^p k v kpg^ a + kpkpg va kpk a g V p iD^{k ) , 



we exploit again the anti-symmetry of both surface elements to obtain 



NP 

A.nc 



tt 2 g 2 (i - k) r . /•'. /• d-k <- r 

3(iV c 2 - 1) 



dc\ I dr-> , .......... 



lim / dui I du 2 e 

T 



-ikz 



2 (n • k) (r 2 • k) - (fx • r 2 ) (k - k s )(k° + k A ) tD[(k 2 ) (2.52) 



with z as given in Q2.46|) . Again the identities % • r 2 = T\ ■ n e = and n@ ■ n Q = 2 
have been used. Performing the u\ and u 2 integrations in the limit T — > 00, one 
obtains — as in Q2.47 ) — two ^-functions which allow us to carry out the integrations 
over k° and k 3 immediately. This leads to 



NP. 



Tl 2 G 2 {l-K) 

3(iV c 2 - 1) 

7T 2 G 2 (1 - K) 

3(iV 2 - 1) 



W fdv 2 (^^^ 
(i Jo J V-k)- 

' 9 (\ 8 -n'( 2 ) 
dvx Jo av 2 



^1 {k!_){n-k ± ){f 2 -k ± )e 



dv lTr - / dv 2Tr - iD[ ( \r lq + v x r x ~ r 2q - v 2 f 2 ) 



ik± (fig +01 fi -fa, -v 2 r 2 ) 



(2.53) 



where iD^ is the non-confining correlation function in transverse space defined 



analogously to (|2.48|) . The v\ and v 2 integrations are trivial and lead (cf. Fig. |]b) 
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to 



NP 
A.nc 



7T 2 Go 



3(N* - 1) 



'1-k) 



iD'} 2) (f lq - r 2q ) + iD[ i2) (fig - r 2q ) 

- lD^ 2) (fig - f 2 q) ~ (fg - f 2 g 



(2.54) 



Using D' x (A;^), derived from the exponential correlation function ( |2.34| ) in Ap- 
pendix B we obtain 



iD[ {2 \z ± ) = 7T a 4 [3 + 3(|4|/a) + (|^|/a) 2 ] exp(-|4| 



^2.55) 



The non-perturbative components, \c P an d X^f > l ea d to co/or transparency for 
small dipoles, i.e. a dipole-dipole cross section with <JDD{ri,r 2 ) oc |fi| 2 |f2| 2 for 
I ^12 1 ~~ > 0, as known for the perturbative case [34]. This can be seen by squar- 



ing ( |2.49| ) and ( j2.54| ) to obtain the leading terms in the T-matrix element for small 
dipoles (see ( p.63| )). 

The perturbative component \ P is defined as 

■2 



X 



Si 



lJa^(x 1 ) I <h'" : (.r 2 )F; wii -,- 2 ; 



s 2 



9_ 
4 



da^ixx) / da pa {x 2 ) 
Si Js 2 



(2.56) 



-ikz 



kuk a 9/j,p k u k p g^ a + k^k p g va k p k a g vp iD p {k ^ 



and shows a structure identical to the one of Xnc given in (|2.52j ) . Accounting for the 
different prefactors and the different correlation function, the result for Xnc ( P-54|) 
can be used to obtain 



X 



g\fi q - r 2q ) iD'^ 2) (r Xq - r 2q ) + £f 2 (fig - f 2q ) %D'^ } {r lq - f 2q ) 



1(2) 



- g\n q - f 2 g) iD'^ ] (fg - f 2 g) ~ g' '' (f^g ~ f 2 g) iD'^' (fig - f 2q 

where the running coupling g 2 (z±) is understood as given in ( |2.39| ). With 
obtains the perturbative correlation function in transverse space 



t(2) 



tD'^ 2) (z ± ) = —K (m G \z ± \) 
Lix 



(2.57) 



one 



(2.58) 



where Kq denotes the modified Bessel function (McDonald function) 
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NP 
nc 



In contrast to the confining component \c F '■> the non-confining components, x 
and \ P i depend only on the transverse position between the quark and antiquark 
of the two dipoles and are therefore independent of the surface choice. 



Finally, we explain that the vanishing of XS1S1 an d XS2S2 anticipated in Sec. |2J 



results from the light-like loops and the tensor structures in F^ up(7 . Concentrat- 
ing — without loss of generality — on XSiSu the appropriate infinitesimal surface 
elements ( |2.43| ) and the F M1/po .-ansatz given in ( |2.31| ), ( |2.32| ), and (|2.37| ) are in- 



serted into ( 2.23 ). Having simplified the resulting expression by exploiting the anti- 



symmetry of the surface elements, one finds only terms proportional to ni 



and n ffi ■ z with z h 



x 1 x 2 



(wi - u 2 )ng + [vi - v 2 )r^ 



Since £ = and 



0, which is evident from (|2.42 ), all terms vanish and XS1S1 = is derived. 



Note that x 
wave functions \ipi{zi, 
the replacement (fj — > — r iy Zi 
imaginary and reads for N c = 3 



Xc P + Xnc + X P is a real-valued function. Since, in addition, the 
fj)| 2 used in this work (cf. Appendix are invariant under 



1 — Z{), the T-matrix element becomes purely 



T(s,t) = 2is Jd 2 b ± e i<LS± Jdz 1 d 2 r 1 Jdz 2 d 2 r 2 



2 (I , r _ ^\ 1 (2 , r 
1 - g cos I -X{b±,Zi,ri,Z2,r 2 ) J - - cos I -x{o±, z u n, z 2 , r 2 



(2.59) 



The real part averages out in the integration over and Zi since the x-function 
changes sign 

X(6±, 1 - zi, -r u z 2} r 2 ) = ~x(b±, z 1 ,r 1 , z 2 ,r 2 ) , (2.60) 



which can be seen directly from fl2.49 ),( 2.54|) and (|2.57|) as (fi — > — fx, z\ — > 1 — Zi) 
implies ri g — > r*ig. In physical terms, (r*j — > — ff , — > 1 — 2^) corresponds to 
charge conjugation i.e. the replacement of each parton with its antiparton and the 
associated reversal of the loop direction. 



Consequently, the T-matrix (|2.59| ) describes only charge conjugation C — +1 
exchange. Since in our quenched approximation purely gluonic interactions are 
modelled, fl2.59|) describes only pomeronP] but not reggeon exchange. 



8 Oddcron C = —1 exchange is excluded in our model. It would survive in the following cases: 

(a) Wave functions are used that are not invariant under the transformation (f\ — > — fi, Zi — * 1 — Zi). 

(b) The proton is described as a system of three quarks with finite separations modelled by three 
loops with one common light-like line, (c) The Gaussian approximation that enforces the truncation 
of the cumulant expansion is relaxed and additional higher cumulants are taken into account. 
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2.4 Energy Dependence 



Until now, the derived T-matrix element leads to energy independent total cross 
sections in contradiction to the experimental observation. In this section, we intro- 
duce the energy dependence in a phenomenological way inspired by other successful 
models. 



Most models for high-energy scattering are constructed to describe either hadron- 
hadron or photon-hadron reactions. For example, Kopeliovich et al. PSJ as well as 



Berger and Nachtmann focus on hadron-hadron scattering. In contrast, Golec- 
Biernat and Wiisthoff || and Forshaw, Kerley, and Shaw |36| concentrate on photon- 
proton reactions. A model that describes the energy dependence in both hadron- 
hadron and photon-hadron reactions up to large photon virtualities is the two- 
pomeron model of Donnachie and Landshoff || . Based on Regge theory, they find a 
soft pomeron trajectory with intercept 1 + e so f t ~ 1.08 that governs the weak energy 
dependence of hadron-hadron or j*p reactions with low Q 2 and a hard pomeron 
trajectory with intercept 1 + e^ard ~ 1.4 that governs the strong energy dependence 
of 7*p reactions with high Q 2 . Similarly, we aim at a simultaneous description 
of hadron-hadron, photon-proton, and photon-photon reactions involving real and 
virtual photons as well. 



In line with other two-component (soft + hard) models || 23|, |24|, |36], |37| and the 
different hadronization mechanisms in soft and hard collisions, our physical ansatz 
demands that the perturbative and non-perturbative contributions do not interfere. 
Therefore, we modify the cosine-summation in ( |2.59| ) allowing only even numbers of 
soft and hard correlations, (x NP ) " (x P ) with n,m G IV. Interference terms with 
odd numbers of soft and hard correlations are subtracted by the replacement 



cos [ex] = cos [c (x NP + X P )] -> cos [c% NP ] cos [cx P ] 



(2.61) 



where c = 1/3 or 2/3. This prescription leads to the following factorization of soft 
and hard physics in the T-matrix element, 



T(s,t) = 2is d 2 b ± e i<kk d Zl d 2 n dz 2 d 2 r 2 |Vi(2i,ri)| 2 \ fo(z 2 , r 2 )\ 



x 



1 ~ I C0 \3 XNP ) C ° S (^ P ) " I C °\3 XNP ) C ° S (^ P 

NP 



(2.62) 



In the limit of small x-functions, x ^ 1 an d X ^ 1> one g e ^s 

T{s,t) = 2is [d 2 b ± e i(L ^ [d Zl d 2 ri [dz 2 d 2 r 2 |Vi(^i,fi)| 2 \ip 2 {z 2 ,f 2 )\ 



(x NP ) 2 +(x P ) 



P\2 



(2.63) 



19 



In this limit, the T-matrix element evidently becomes a sum of a perturbative and 
a non-perturbative component. Of course, the perturbative component, (x P ) 2 ; co- 
incides with the well-known perturbative two-gluon exchange [ 17| . Correspondingly, 



the non-perturbative component, (x NP ) 2 , represents the non-perturbative gluonic 
interaction on the "two-gluon-exchange" level. 

As the two-component structure of ( |2.63| ) reminds of the two-pomeron model of 
Donnachie and Landshoff H], we adopt the powerlike energy increase and ascribe a 
weak energy dependence to the non-perturbative component \ NP an d a strong one 
to the perturbative component \ P 

WP 

(X NP ) 2 - {x NP {s)) 2 :={x NP ) 2 {-^ 

p 

with the scaling factor sqRq. The power like energy dependence with the exponents 
e NP < e p < 1 guarantees Regge type behavior at moderately high energies, 
where the small- x limit (|2.63|) is appropriate. In (|2.64|) , the energy variable s is 



scaled by the factor r-j 2 that ahows to rewrite the energy dependence in photon- 
hadron scattering in terms of the appropriate Bjorken scaling variable x 



,f 1 2 oc^ = ^, (2.65) 
where |fi| is the transverse extension of the qq dipole in the photon. A similar 



factor has been used before in the dipole model of Forshaw, Kerley, and Shaw |3q] 



and also in the model of Donnachie and Dosch |37] in order to respect the scaling 
properties observed in the structure function of the proton.^] In the dipole-proton 
cross section of Golec-Biernat and Wiisthoff ||, Bjorken x is used directly as energy 
variable which is important for the success of the model. In fact, also in our model, 
the f*i ff factor improves the description of 7*p reactions at large Q 2 . 

The powerlike Regge type energy dependence introduced in (|2.64j) is, of course, 
not mandatory but allows successful fits and can also be derived in other theoreti- 
cal frameworks: A powerlike energy dependence is found for hadronic reactions by 
Kopeliovich et al. ]35| and for hard photon-proton reactions from the BFKL equa- 
tion M. However, these approaches need unitarization since their powerlike energy 



9 In the model of Donnachie and Dosch f37|| , s |fi| |f%| is used as the energy variable if both 
dipoles are small, which is in accordance with the choice of the typical BFKL energy scale but 
leads to discontinuities in the dipolc-dipole cross section. In order to avoid such discontinuities, 



we use the energy variable ( 2.64 ) also for the scattering of two small dipoles. 
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dependence will ultimately violate S-matrix unitarity at asymptotic energies. In our 
model, we use the following T-matrix element as the basis for the rest of this work 



T{s,t) = 2is d 2 b ± e 1 ^ \dz x d 2 r x dz 2 d 2 r 2 l^ifa, ?i)| 2 \M^,r 2 )\ 



1 cos 

3 



X^(s) cos ^(s) 



cos 



X Ny (s))cos(^-x P (s) 



[2.66) 



where the cosine functions ensure the unitarity condition in impact parameter space 
as shown in Sec. Indeed, the multiple gluonic interactions associated with the 
higher order terms in the expansion of the cosine functions are important for the 
saturation effects observed within our model at ultra-high energies. 

Having ascribed the energy dependence to the x-function, the energy behavior 
of hadron-hadron, photon-hadron, and photon-photon scattering results exclusively 
from the universal loop-loop correlation function Sdd- 



2.5 Model Parameters 

Lattice QCD simulations provide important information and constraints on the 
model parameters. The fine tuning of the parameters was, however, directly per- 
formed on the high-energy scattering data for hadron-hadron, photon-hadron, and 
photon-photon reactions where an error (x 2 ) minimization was not feasible because 
of the non-trivial multi-dimensional integrals in the T-matrix element fl2.66|) . 

The parameters a, k, G 2 , ma, M 2 , soRq, e and e p determine the dipole- 
dipole scattering and are universal for all reactions described. In addition, there 
are reaction-dependent parameters associated with the wave functions which are 
provided in Appendix 

The non-perturbative component involves the correlation length a, the gluon 
condensate G 2 , and the parameter k indicating the non-Abelian character of the 
correlator. With the simple exponential correlation functions specified in Euclidean 
space-time ( |2.34j ), we obtain the following values for the parameters of the non- 
perturbative correlator (|2.32|) 



a = 0.302 fm, k = 0.74, G 2 = 0.074 GeV 4 , (2.67) 
and, correspondingly, the string tension 

a = ^ K ° 2 ^ = 0.22 GeV 2 = 1.12 GeV/fm , (2.68) 
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which is consistent with hadron spectroscopy f55fl , Regge theory P5[ , and lattice 
QCD investigations f3l |. 



Lattice QCD computations of the gluon field strength correlator down to dis- 
tances of 0.4 fm have obtained the following values with the exponential correlation 
function (fZ3|) |TJ: a = 0.219 fm, k = 0.746, G 2 = 0.173 GeV 4 . This value for k 



is in agreement with the one in ( 2.67 ), while the fit to high-energy scattering data 
clearly requires a larger value for a and a smaller value for G^- 

The perturbative component involves the gluon mass vtiq as IR regulator (or 
inverse "perturbative correlation length") and the parameter M 2 that freezes the 
running coupling ( |2.39|) for large distance scales at the value a s = 0.4, where the 
non-perturbative component of our model with the above ingredients is at work 



according to a low energy theorem |16| , |33| . We adopt the parameters 
m G = m p = 0.77 GeV and M 2 = 1.04GeV 2 . 



(2.69) 



The energy dependence of the model is associated with the energy exponents 
e NP and e p , and the scaling parameter s Rq 

6^ = 0.125, e p = 0.73, and s R* = (47GeVfm 2 ) 2 . (2.70) 

In comparison with the energy exponents of Donnachie and Landshoff ||[7]], e so ft ~ 
0.08 and €h ar d ~ 0.4, our exponents are larger. However, the cosine functions in our 
T-matrix element ( |2.66| ) reduce the large exponents so that the energy dependence 
of the cross sections agrees with the experimental data as illustrated in Sec |. 



3 Impact Parameter Profiles and S'-Matrix Uni- 
tarity 

In this section, the S'-matrix unitarity is analysed in our model. On the basis of the 
impact parameter dependence of the scattering amplitude, saturation effects can be 
exposed that manifest the unitarity of the S'-matrix. For each impact parameter 
the energy at which the unitarity limit becomes important can be determined. This 
is used to show the saturation of the gluon distribution and to localize saturation 
effects in experimental observables. 

The impact parameter dependence of the scattering amplitude is given by T(s, 

T(s,t = -q 2 ) = 4s fd 2 b ± e i ^T(s,\b ± \) (3.1) 
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and in particular by the profile function 

J(s,\b ± \) = 2lmT(s,\b ± \) , (3.2) 

which describes the blackness or opacity of the interacting particles as a function of 
the impact parameter and the cm. energy ^/s. In fact, the profile function ( |3.2|) 
determines every observable if the T-matrix is — as in our model — purely imagi- 
nary. 

The S'-matrix unitarity, SS' = S^S = 11, leads directly to the unitarity condition 
in impact parameter space ] 

ImT(s, |6l|) = \T(s, \b ± \)\ 2 + G inel ( Sj \b ± \) , (3.3) 

where G ine i(s, \b±\) > is the inelastic overlap function ^2|.Q This unitarity condi- 
tion imposes an absolute limit on the profile function 

< 2\T(s,\b ± \)\ 2 < J(s,\b ± \) < 2 (3.4) 

and the inelastic overlap function, G ine i(s, \b±\) < 1/4. At high energies, however, 
the elastic amplitude is expected to be purely imaginary. Consequently, the solution 
of (|3.3|) reads 



J(s,\b ± \) = l±^l-4G ine »(s,|&]_|) (3.5) 

and leads with the minus sign corresponding to the physical situation to the reduced 
unitarity bound 

< J(s, \b±\) < 1 . (3.6) 

Reaching the black disc limit or maximum opacity at a certain impact parameter 
\b±\, J(s, \b±\) = 1, corresponds to maximal inelastic absorption G ine i(s, \b±\) = 1/4 
and equal elastic and inelastic contributions to the total cross section at that impact 
parameter. 

In our model, every reaction is reduced to dipole-dipole scattering with well 
defined dipole sizes \f\\ and longitudinal quark momentum fractions Zi. The unitarity 
condition in our model becomes, therefore, most explicit in the profile function 

Jdd(s, \b±\,zi, \ri\,z 2 , |f 2 |) = / ^p- [ ^p- 1 - S DD (s, b±, z u f u z 2 ,r 2 ) , (3.7) 

J 2ir J 2ir L J 

where fa describes the dipole orientation, i.e. the angle between f\ and b±, and Sdd 
describes elastic dipole-dipole scattering 

Sdd = I^Qx^to) cosQx P (^)) + \ cosQx^w) cosQ X P ( S )) (3.8) 



10 Integrating (3.3) over the impact parameter space and multiplying by a factor of 4 one obtains 
the relation a tot {s) = a el {s) + a mel (s). 
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with the purely real- valued eikonal functions x NP { s ) an d X P { S ) defined in ( |2.64| ). 
Because of \Sdd\ < 1, a consequence of the cosine functions in (|3.8|) describing 
multiple gluonic interactions, J^d respects the absolute limit (|3~4l ). Thus, the elastic 



dipole-dipole scattering respects the unitarity condition ( |3.3| ). At high energies, the 
arguments of the cosine functions in Sod become so large that these cosines average 
to zero in the integration over the dipole orientations. This leads to the black disc 
limit Jpff = 1 reached at high energies first for small impact parameters. 

If one considers the scattering of two dipoles with fixed orientation, the inelastic 
overlap function obtained from the unitarity constraint 

^DD 

/ -, 2 / 1 /vp /x\ /l p / \ \ 1 /2 , N \ / 2 p , 

S 




jjcosQx^( S )) cosQx P (s)) + IcosQx^Oo) cosQx P ( 




shows nonphysical behavior with increasing energy. This behavior is a consequence 
of aritifically fixing the orientations of the dipoles. If one averages over the dipole 
orientations as in all high-energy reactions considered in this work, no unphysical 
behavior is observed. 



3.1 The Profile Function for Proton-Proton Scattering 

The profile function for proton-proton scattering 

1 - S DD (s,b±,z 1 ,r 1 ,z 2 ,r 2 



J PP {s, \b±\) = Jdzid 2 ri J dz 2 d 2 r 2 \i! p (zi,fi)\ 2 \ifj p (z 2 ,f 2 )\ 2 

(3.10) 

is obtained from (|3.7|) by weighting the dipole sizes \fi\ and longitudinal quark 
momentum fractions Z{ with the proton wave function \ifj p (zi, fi)\ 2 from Appendix [A]. 



Using the model parameters from Sec. one obtains the profile function 
Jpp( s , \b±\) shown in Fig. ^| for cm. energies from y/s = 10 GeV to y/s = 10 8 GeV. 
Up to \fs ~ 100 GeV, the profile has approximately a Gaussian shape. Above 
y/s = 1 TeV, it significantly develops into a broader and higher profile until the 
black disc limit is reached for y/s ~ 10 6 GeV and |6j_| =0. At this point, the cosine 
functions in Sdd average to zero 

JdzKPn Jdz 2 d 2 r 2 \^ p (z 1 ,n)\ 2 \M^,r 2 )\ 2 S DD (^> 10 6 GeV, \b L \ = 0, . . .) w 

(3.11) 
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Figure 2: The profile function for proton-proton scattering J pp (s, \b±\) is shown versus the 
impact parameter \b±\ for cm. energies from y/s = 10 GeV to y/s = 10 8 GeV. The unitarity 



limit (3.4) corresponds to J pp (s, \b±\) = 2 and the black disc limit (pTq) to J pp (s, \b±\) = 1. 



so that the proton wave function normalization determines the maximum opacity 

J™ x = Jd Zl d 2 ri Jdz 2 d 2 r 2 l^i.fOI 2 |^2,r 2 )| 2 = 1 . (3.12) 

Once the maximum opacity is reached at a certain impact parameter, the profile 
function saturates at that and extends towards larger impact parameters with 
increasing energy. Thus, the multiple gluonic interactions important to respect the 
S'-matrix unitarity constraint (|3.3j) lead to saturation for y/s > 10 6 GeV. 



The above behavior of the profile function illustrates the evolution of the proton 
with increasing cm. energy. The proton is gray and of small transverse size at 
small y/s but becomes blacker and more transversally extended with increasing y/s 
until it reaches the black disc limit in its center at y/s w 10 6 GeV. Beyond this 
energy, the proton cannot become blacker in its central region but in its periphery 
with continuing transverse growth. Furthermore, the proton boundary seems to stay 



diffusive as claimed also in 13 . 
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According to our model the black disc limit will not be reached at LHC. Our 
prediction of y/s ~ 10 6 GeV = 10 3 TeV for the onset of the black disc limit in 
proton-proton collisions is about two orders of magnitude beyond the LHC energy 
y/s = 14TeV. This is in contrast, for example, with where the value predicted 



for the onset of the black disc limit is ^fs = 2 TeV, i.e. small enough to be reached at 
LHC. However, we feel confidence in our LHC prediction since our profile function 
Jpp(s, \b±\) yields good agreement with experimental data for cross sections up to 
the highest energies as shown in Sec. [5| 

For hadron-hadron reactions in general, the wave function normalization of the 
hadrons determines the maximum opacity analogous to ( |3.12|) and the transverse 
hadron size the cm. energy at which it is reached. Consequently, the maximum 
opacity obtained for np and Kp scattering is identical to the one for pp scattering 
due to the normalization (|A.2| ). Furthermore, the smaller size of pions and kaons 
in comparison to protons demands slightly higher cm. energies to reach this maxi- 
mum opacity. This size effect becomes more apparent in longitudinal photon-proton 
scattering, where the size of the dipole emerging from the photon can be controlled 
by the photon virtuality. 



3.2 The Profile Function for Photon-Proton Scattering 

The profile function for a longitudinal photon 7^ scattering off a proton p 

Jr L p( s Ab±\,Q 2 ) = [dz 1 d 2 r 1 Jdz 2 d 2 r 2 \^* L (z 1 ,r 1 ,Q 2 )\ 2 \4> p (z 2 ,r 2 )\ 2 

1 - S DD (b±,s, zi,fi,z 2 ,f 2 ) (3.13) 

is calculated with the longitudinal photon wave function \i/) 7 *(zi, fi, Q 2 )\ 2 given 
in ( |A.5| ). In this way, the profile function ( |3.13|) is ideally suited for the investi- 



gation of dipole size effects since the photon virtuality Q 2 determines the transverse 
size of the dipole into which the photon fluctuates before it interacts with the proton. 

Figure [| shows the \b±\ dependence of the profile function J^ p (s, \b±\, Q 2 ) divided 
by a/ir for cm. energies -^i from 10 GeV to 10 9 GeV and a photon virtuality of Q 2 = 
1 GeV 2 , where a is the fine-structure constant. One clearly sees that the qualitative 
behavior of this rescaled profile function is similar to the one for proton-proton 
scattering. However, the black disc limit induced by the underlying dipole-dipole 
scattering depends on the photon virtuality Q 2 and is given by the normalization of 
the longitudinal photon wave function 

J™p(Q 2 ) = Jdzd 2 r\^ L (z,r,Q 2 )\ 2 (3.14) 
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Figure 3: The profile function for a longitudinal photon scattering off a proton 
^7lp( s ' I^J-li Q 2 ) divided by ajix is shown versus the impact parameter |6j_| at a photon 
virtuality of Q 2 = 1 GeV 2 and cm. energies from y/s = 10 GeV to y/s = 10 9 GeV. The 



value of the black disc limit is J™p X (Q 2 



lGeV z 



0.00164 



since the proton wave function is normalized to one. 

The photon virtuality Q 2 does not only determine the absolute value of the 
black disc limit but also the cm. energy at which it is reached. This is illustrated in 
Fig. |, where the y/s dependence of J 7 * p (s, \b±\ = 0, Q 2 ) divided by ch/k is presented 
for Q 2 = 1, 10, and 100 GeV 2 . With increasing resolution Q 2 , i.e. decreasing dipole 
sizes, \r 7 ^\ 2 oc l/Q 2 , the absolute value of the black disc limit grows and higher 
energies are needed to reach this limit .0 The growth of the absolute value of the 
black disc limit is simply due to the normalization of the longitudinal photon wave 
function while the requirement of higher energies to reach this limit is due to the 
decreasing interaction strength with decreasing dipole size. The latter explains also 
why the energies needed to reach the black disc limit in and Kp scattering are 



11 Note that the Bjorken x at which the black disc limit is reached decreases with increasing 
photon virtuality Q 2 . (See also Fig. ||) 
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Figure 4: The profile function for a longitudinal photon scattering off a proton 
J^* L p{ s -> I&lIj Q 2 ) divided by a/ir is shown versus the cm. energy ^fs at zero impact param- 
eter (|fe|_| = 0) for photon virtualities Q 2 = 1, 10, and 100 GeV 2 . 



higher than in pp scattering. Comparing ^* L p scattering at Q 2 = 1 GeV 2 with pp 
scattering quantitatively, the black disc limit J™ a p x (Q 2 = 1 GeV 2 ) = 0.00164 is about 
three orders of magnitude smaller because of the photon wave function normalization 
(oc a/n). At \b±\ = it is reached at an energy of ^fs ~ 10 8 GeV, which is about 
two orders of magnitude higher because of the smaller dipoles involved. 

The way in which the profile function J 7 £ P (s, |6j_|, Q 2 ) approaches the black disc 
limit at high energies depends on the shape of the proton and longitudinal photon 
wave function at small dipole sizes |r 12 |. At high energies, dipoles of typical sizes 
< |ri j2 | < i?o ( s o/ s ) 1//4 give the main contribution to S^ p = 1 — J 7 £ P because 
of ( p. 64 ) and the fact that the contribution of the large dipole sizes averages to zero 
upon integration over the dipole orientations, cf. also (|3.11|) . Since S y * p is a measure 
of the proton transmittance, this means that only small dipoles can penetrate the 
proton at high energies. Increasing the energy further, even these small dipoles are 
absorbed and the black disc limit is reached. However, the dependence of the profile 
function on the short distance behavior of normalizable wave functions is weak which 
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can be understood as follows. Because of color transparency, the eikonal functions 
X NP (s) and x P ( s ) are small for small dipole sizes < |r li2 | < Ro (so/s) 1 ^ at large 
y/s. Consequently, Sdd ~ 1 and 

Jr L p( s i \b±\,Q 2 ) 

1 r c (s) 1 r c (s) 

JT*(Q 2 )-47T 2 [dzi [dr 1 r 1 \^(z 1 ,r 1 ,Q 2 )\ 2 [dz 2 [dr 2 r 2 \ij p (z 2 ,r 2 )\ 2 (3.15) 





where r c (s) ~ i?o (so/s) 1 ^- Clearly, the linear behavior from the phase space factors 
ri 5 2 dominates over the r\ ^-dependence of normalizable wave functions] 17 ! More 
generally, for any profile function involving normalizable wave functions, the way in 
which the black disc limit is approached depends only weakly on the short distance 
behavior of the wave functions. 



4 A Scenario for Gluon Saturation 



In this section, we estimate the impact parameter dependent gluon distribution of 
the proton xG(x,Q 2 , \b±\). Using a leading twist, next-to-leading order QCD rela- 
tion between xG(x, Q 2 ) and the longitudinal structure function F L (x, Q 2 ), we relate 
xG(x,Q 2 , \b±\) to the profile function J 7 t p (s = Q 2 /x, \b±\,Q 2 ) and find low-x sat- 
uration of xG(x,Q 2 , \b±\) as a manifestation of S'-matrix unitarity. The resulting 
xG(x, Q 2 , \b±\) is, of course, only an estimate since our profile function contains also 
higher twist contributions. Furthermore, in the considered low-x region, the lead- 
ing twist, next-to-leading order QCD formula may be inadequate as higher twist 



contributions |45[ and higher order QCD corrections |4q,|47| are expected to be- 
come important. Nevertheless, still assuming a close relation between Fl(x,Q 2 ) 
and xG(x, Q 2 ) at low x, we think that our approach provides some insight into the 
gluon distribution as a function of the impact parameter and into its saturation. 

The gluon distribution of the proton xG(x, Q 2 ) has the following meaning: 
xG(x, Q 2 )dx gives the momentum fraction of the proton which is carried by the glu- 
ons in the interval [x, x+dx] as seen by probes of virtuality Q 2 . The impact parameter 
dependent gluon distribution xG(x, Q 2 , \b±\) is the gluon distribution xG(x, Q 2 ) at a 



12l 



For our choice of the wave functions in (3.15), one sees very explicitly that the specific Gaus- 



sian behavior of \tp p (z2, ^2)| 2 and the logarithmic short distance behavior of \ip-yt (zi, ri, Q 2 )\ 2 is 
dominated by the phase space factors ri : 2- 
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given impact parameter so that 

In leading twist, next-to-leading order QCD, the gluon distribution xG(x, Q 2 ) is 
related to the structure functions Fl(x,Q 2 ) and F 2 (x, Q 2 ) of the proton 



F L {x,Q 2 



71 



4 ( x dy ( x 
3J X y 



(4-2) 

where Y2f e / is a flavor sum over the quark charges squared. For four active flavors 



and x < 10 , ( |4.2| ) can be approximated as follows fl9| 



xG(x,Q 2 ) « -5.8 
5 



' U F L (0A17x,Q 2 ) --]-F 2 (0.75x,Q 2 ) 



4a, " v 7 1.97 



(4.3) 



For typical A QCD = 100 - 300 MeV and Q 2 = 50 - 100 GeV 2 , the coefficient of 
Fl in (|4.3j), 37r/(4a s ) = O(10), is large compared to the one of F 2 . Taking into 
account also the values of F 2 and F L , in this Q 2 region and for x < 10~ 3 , the gluon 
distribution is mainly determined by the longitudinal structure function. The latter 
can be expressed in terms of the profile function for longitudinal photon-proton 
scattering using the optical theorem (cf. (|5.1|) ) 



where the s-dependence of the profile function is rewritten in terms of the Bjorken 
scaling variable, x = Q 2 /s. Neglecting the F 2 term in (f4.3p, consequently, the gluon 
distribution reduces to 

xG(x, Q 2 ) « 1.305 - fd 2 b ± J 7 «„(0.417x, \b±\,Q 2 ) . (4.5) 
n 2 a s a J L 



Comparing (|4.1| ) with (|4.5|) , it seems natural to relate the integrand of (|4.5|) to the 
impact parameter dependent gluon distribution 

O 2 IT 

xG(x, Q 2 , \b ± \) « 1.305 -J J 7 *„(0.417x, \b ± \,Q 2 ) . (4.6) 

7T (x s a L 

The black disc limit of the profile function for longitudinal photon-proton scat- 
tering ( |3.14j ) imposes accordingly an upper bound on xG(x,Q 2 , \b±\) 



xG(x,Q 2 ,\b ± \) < xG max (Q 2 ) « 1.305 -f- - JT P X (Q 2 ) , (4.7) 
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which is the low- a; saturation value of the gluon distribution xG(x,Q 2 , \b±\) in our 
approach. With tij 
of (14.71) is obtained 



approach. With ir J™ a p (Q 2 ) / a ~ 1 as shown in Fig. f|, a compact approximation 



xG(x,Q 2 , \b ± \) < xG max (Q 2 ) « , (4.8) 



IX Ot 



which is consistent with the results in [|47|, p0|, plf and indicates strong color-field 
strengths G a ~ 1/ \fcT s as well. 



According to our relations (^4.6|) and (^4.7|) , the blackness described by the pro- 
file function is a measure for the gluon distribution and the black disc limit corre- 
sponds to the maximum gluon distribution reached at the impact parameter under 
consideration. In accordance with the behavior of the profile function J-y* p , see 
Fig. [I the gluon distribution xG(x,Q 2 , \b±\) decreases with increasing impact pa- 
rameter for given values of x and Q 2 . The gluon density, consequently, has its 
maximum in the geometrical center of the proton, i.e. at zero impact parameter, 
and decreases towards the periphery. With decreasing x at given Q 2 , the gluon 
distribution xG(x, Q 2 , \b±\) increases and extends towards larger impact parameters 
just as the profile function J 7 £ P for increasing s. The saturation of the gluon dis- 
tribution xG(x,Q 2 , \b±\) sets in first in the center of the proton (\b±\ = 0) at very 
small Bjorken x. 

In Fig. |5], the gluon distribution xG(x, Q 2 , \b±\ = 0) is shown as a function 
of x for Q 2 = 1, 10, andlOOGeV 2 , where the relation (|4.6|) has been used also 
for low photon virtualities. Evidently, the gluon distribution xG(x, Q 2 , \b±\ = 0) 
saturates at very low values of x < 10~ 10 for Q 2 > 1 GeV 2 . The photon virtuality 
Q 2 determines the saturation value ( |4.7|) and the Bjorken-x at which it is reached 
(cf. also Fig. [3]). For larger Q 2 , the low-x saturation value is larger and is reached at 
smaller values of x, as claimed also in [Q. Moreover, the growth of xG(x, Q 2 , \b±\ = 
0) with decreasing x becomes stronger with increasing Q 2 . This results from the 
stronger energy increase of the perturbative component, e p = 0.73, that becomes 
more important with decreasing dipole size. 

According to our approach, the onset of the xG(x,Q 2 , |&j_|)-saturation appears 
for Q 2 > 1 GeV 2 at x < 10 -10 , which is far below the x-region accessible at HERA 
(x > 10~ 6 ). Even for THERA (x > 10~ 7 ), gluon saturation is not predicted for 
Q 2 > 1 GeV 2 . However, since the HERA data can be described by models with and 
without saturation embedded f52|, the present situation is not conclusive.^ 



13 So far, the most striking hint for saturation in the present HERA data at x ~ 10 4 and 
Q 2 < 2 GeV 2 has been the turnover of dF2(x, Q 2 )/dln(Q 2 ) towards small x in the Caldwell plot [p3| , 
which is still a controversial issue due to the correlation of Q 2 and x values. 
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Figure 5: The gluon distribution of the proton at zero impact parameter xG(x, Q 2 , \b±\ = 
0) is shown as a function of x for Q 2 = 1, 10, andlOOGeV 2 . The results are obtained 
within the approximation fl4.6|) . 

Note that the S'-matrix unitarity condition ( |3.3| ) together with ( |4.6|) requires the 
saturation of the impact parameter dependent gluon distribution xG(x, Q 2 , \b±\) 
but not the saturation of the integrated gluon distribution xG(x, Q 2 ). Due to 
multiple gluonic interactions in our model, this requirement is fulfilled, as can be 
seen from Fig. |3| and relation ( |4.6|) . Indeed, approximating the gluon distribution 
xG(x, Q 2 , \b±\) in the saturation regime of very low a; by a step-function 

xG(x,Q 2 ,\b ± \)^xG max (Q 2 )Q(R(x,Q 2 )-\b ± \) , (4.9) 

where R(x, Q 2 ) denotes the full width at half maximum of the profile function, one 
obtains with (|4.1| ), ( [17D and (|4.8| ) the integrated gluon distribution 

xG(x ^) „ 1M5 Q 2R2 ^l J?7m „ Q 2R2 M 2 ) , (4.io) 

7ia s a 7lP 7ia s 

which does not saturate because of the increase of the effective proton radius R(x, Q 2 ) 
with decreasing x. Nevertheless, although xG(x, Q 2 ) does not saturate, the satura- 
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tion of xG(x,Q 2 , \b±\) leads to a slow-down in its growth towards small x.f^\ Inter- 
estingly, our result ( |4. 10| ) coincides with the result of Mueller and Qiu 



Finally, it must be emphasized that the low-x saturation of xG(x,Q 2 , \b±\), re- 
quired in our approach by the S'-matrix unitarity, is realized by multiple gluonic in- 
teractions. In other approaches that describe the evolution of the gluon distribution 
with varying x and Q 2 , gluon recombination leads to gluon saturation [f46| , f|7], [54H56 



which is reached when the probability of a gluon splitting up into two is equal to the 
probability of two gluons fusing into one. A more phenomenological understanding 
of saturation is attempted in P, |5T . 



5 Comparison with Experimental Data 

In this section, we discuss the phenomenological performance of our model. We com- 
pute total, differential, and elastic cross sections, structure functions, and diffractive 
slopes for hadron-hadron, photon-proton, and photon-photon scattering, compare 
the results with experimental data including cosmic ray data, and provide predic- 
tions for future experiments. Having studied the saturation of the impact parameter 
profiles, we show here how this manifestation of unitarity translates into the quan- 
tities mentioned above and how it could become observable. 



Using the T-matrix ( |2.66|) with the parameters and wave functions from Sec. |2] 



and Appendix [X], we compute the pomeron contribution to pp, pp, n^p, K ± p, 7*p, 
and 77 reactions in terms of the universal dipole-dipole scattering amplitude Sod- 
This allows one to compare reactions induced by hadrons and photons in a systematic 
way. In fact, it is our aim to provide a unified description of all these reactions and 
to show in this way that the pomeron contribution to the above reactions is universal 
and can be traced back to the dipole-dipole scattering amplitude Sdd- 

Our model describes pomeron (C = +1 gluon exchange) but neither odderon 
(C = — 1 gluon exchange) nor reggeon exchange (quark-antiquark exchange) as 



discussed in Sec. |2.3[ Only in the computation of the hadronic total cross sections 
the reggeon contribution is added ]7], [58] . This improves the agreement with the data 
for a/s < 100 GeV and describes exactly the differences between ab and ab reactions. 

The fine tuning of the model and wave function parameters was performed on the 
data shown below. The resulting parameter set given in Sec. |275| and Appendix |A| 



14 This is analogous to the rise of the total pp cross section with growing cm. energy that slows 
down as the corresponding profile function J pp (s, \bj_\) reaches its black disc limit as shown in 
Sec. O 
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is used throughout this paper. 



5.1 Total Cross Sections 



The total cross section for the high-energy reaction ab — > X is related via the optical 
theorem to the imaginary part of the forward elastic scattering amplitude and can 
also be expressed in terms of the profile function 



<7*J*(s) = ^ImT(M = 0) = 2 Jd 2 b ± J ab (s,\b ± \) , (5.1) 

where a and b label the initial particles whose masses were neglected as they are 
small in comparison to the cm. energy ^fs. 

We compute the pomeron contribution to the total cross section, cr^ ,IP {s), from 



the T-matrix Q2.66|) , as explained above, and add only here a reggeon contribution 



of the form M, Eg] 

, , m / s \ -0-4525 

^' R (s) = X ab (y^s) , (5.2) 

where X ab depends on the reaction considered: X pp = 56.08 mb, X p p = 98.39 mb, 
X^+p = 27.56 mb, X n - p = 36.02 mb, X K + P = 8.15 mb, X K - p = 26.36 mb, X lp = 
0.129 mb, and X 77 = 605 nb. Accordingly, we obtain the total cross section 

o%(s)=o2> F (8) + o%*(8) (5.3) 
for pp, pp, n ± p, K ± p, 7p and 77 scattering. 

The good agreement of the computed total cross sections with the experimental 
data is shown in Fig. |6]. Here, the solid lines represent the theoretical results for 
pp, Tr + p, K + p, 7p, and 77 scattering and the dashed lines the ones for pp, n~p, and 
K~p scattering. The pp, pp, 7r ± p, K ± p, jp |j] and 77 data |59[ taken at accelerators 
are indicated by the closed circles while the closed squares (Fly's eye data) |60H 
and the open circles (Akeno data) |6]J indicate cosmic ray data. Concerning the 



photon-induced reactions, only real photons are considered which are, of course, 
only transverse polarized. 



The prediction for the total pp cross section at LHC (y^i = 14TeV) is a, 



tot 

pp 



114.2 mb in good agreement with the cosmic ray data. Compared with other works, 
our LHC prediction is close to the one of Block et al. UB"2] , cr*°* = 108 ± 3.4 mb, but 



considerably larger than the one of Donnachie and Landshoff 0, o pp = 101.5 mb. 



The differences between ab and ab reactions for \fs < 100 GeV result solely from 
the different reggeon contributions which die out rapidly as the energy increases. The 
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Figure 6: The total cross section a tot is shown as a function of the cm. energy yfs 
for pp, pp, 7r ± p, K^p, jp and 77 scattering. The solid lines represent the model results 
for pp, tt + p, K + p, 7p and 77 scattering and the dashed lines the ones for pp, ir~p, and 
K~p scattering. The pp, pp, tt^p, K^p, jp and 77 data |5{J taken at accelerators are 
indicated by the closed circles while the closed squares (Fly's eye data) |6(| and the open 
circles (Akeno data) indicate cosmic ray data. The arrows at the top point to the 
LHC energy, y/s = 14TeV, and to the onset of the black disc limit in pp (pp) reactions, 
y/s » 10 6 GeV. 
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pomeron contribution to ab and ab reactions is, in contrast, identical and increases 
as the energy increases. It thus governs the total cross sections for yfs > 100 GeV 
where the results for ab and ab reactions coincide. 

The differences between pp (pp) , 7r ± p, and K p scattering result from the differ- 
ent transverse extension parameters, S p = 0.86 fm > S n = 0.607 fm > Sk = 0.55 fm, 
cf. Appendix |A} Since a smaller transverse extension parameter favors smaller 
dipoles, the total cross section becomes smaller, and the short distance physics 
described by the perturbative component becomes more important and leads to a 
stronger energy growth due to e p = 0.73 > e = 0.125. In fact, the ratios cr*°' /er£°* 
and (Tppf&Kp converge slowly towards unity with increasing energy as can already 
be seen in Fig. |[ 



For real photons, the transverse size is governed by the constituent quark masses 
rrif(Q 2 = 0), cf. Appendix [A], where the light quarks have the strongest effect, i.e. 
a*"* oc l/m 2 ld and <r*°' oc 1/m^. Furthermore, in comparison with hadron-hadron 
scattering, there is the additional suppression factor of a for 7p and a 2 for 77 
scattering coming from the photon-dipole transition. In the 77 reaction, also the 
box diagram contributes [^8|,^3| but is neglected since its contribution to the total 
cross section is less than 1% |3^ . 



It is worthwhile mentioning that total cross sections for pp (pp), 7r ± j>, and K ± p 
scattering do not depend on the longitudinal quark momentum distribution in the 
hadrons since the underlying dipole-dipole cross section is independent of the lon- 
gitudinal quark momentum fraction Z{ for t = 0. We show this analytically on the 
two-gluon-exchange level in ||17|| . 

Saturation effects as a manifestation of the S-matrix unitarity can be seen in 
Fig. []. Having investigated the profile function for pp (pp) scattering, we know 
that this profile function becomes higher and broader with increasing energy until 
it saturates the black disc limit first for zero impact parameter (\b±\ = 0) at y/s ~ 
10 6 GeV. Beyond this energy, the profile function cannot become higher but expands 
towards larger values of |6j_|. Consequently, the total cross section (|5.1| ) increases 
no longer due to the growing blackness at the center but only due to the transverse 
expansion of the hadrons. This tames the growth of the total cross section as can 
be seen for the total pp cross section beyond a/s ~ 10 6 GeV in Fig. |6|. 

At energies beyond the onset of the black disc limit at zero impact parameter, 
the profile function can be approximated by 

JT° X {s, |6l|) =N a N b e (r(s) - |6l|) (5.4) 
where N a b denotes the normalization of the wave functions of the scattered particles 
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and R(s) the full width at half maximum of the exact profile function J a b{s, \b±\) that 
reflects the effective radii of the interacting particles. Thus, the energy dependence 
of the total cross section ( f5.ll) is driven exclusively by the increase of the transverse 
extension of the particles R(s) 

o%(s) = 27tN a N b R(s) 2 , (5.5) 



which is known as geometrical scaling The growth of R(s) can at most 

be logarithmic for y/s — > oo because of the Froissart-Lukaszuk-Martin bound ||. 
In fact, a transition from a power-like to an ln 2 -increase of the total pp cross sec- 
tion seems to set in at about y/s ~ 10 6 GeV as visible in Fig. Moreover, since 
the hadronic cross sections join for y/s — > oo, R(s) becomes independent of the 
hadron-hadron reaction considered at asymptotic energies as long as N^f, = 1. Also 
for photons of different virtuality Q\ and Q\ one can check that the ratio of the 
total cross sections ^^(Q 2 ) / a^fJQl) converges to unity at asymptotic energies in 
agreement with the conclusion in |56| . 



5.2 The Proton Structure Function 

The total cross section for the scattering of a transverse (T) and longitudinally 
(L) polarized photon off the proton, <T*|^ p (x, Q 2 ), at photon virtuality Q 2 and cm. 

energy^ squared, s = Q 2 /x, is equivalent to the structure functions of the proton 

F T , L (x,Q 2 ) = ^^f LP (x,Q 2 ) (5.6) 

and 

F 2 (x, Q 2 ) = F T (x, Q 2 ) + F L (x, Q 2 ) . (5.7) 

Reactions induced by virtual photons are particularly interesting because the 
transverse separation of the quark-antiquark pair that emerges from the virtual 
photon decreases as the photon virtuality increases (cf. Appendix |A|) 

2 

\r 7 \ , (5.8) 



Q 2 + Am 2 



u.d 



where m u ^ is a mass of the order of the constituent w-quark mass. With increasing 
virtuality, one probes therefore smaller transverse distance scales of the proton. 

15 Here, \fs refers to the cm. energy in the j*p system. 
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In Fig. |7|, the Q 2 -dependence of the total j*P cross section 

°$ p {s, Q 2 ) = o%,(8, Q 2 ) + o$Ja, Q 2 ) (5.9) 

is presented, where the model results (solid lines) are compared with the experi- 
mental data for cm. energies from yfs = 20 GeV up to ^fs = 245 GeV. Note the 
indicated scaling factors at different yfs values. The low energy data at y/s = 20 GeV 
are from [^7] while the data at higher energies have been measured at HERA by the 
HI |H and ZEUS collaboration §9|. At Q 2 = 0.012 GeV 2 , also the photoproduction 



(Q 2 = 0) data from |7(| are displayed. 



In the window shown in Fig. [7|, the model results are in reasonable agreement with 
the experimental data. The total 7*p cross section levels off towards small values 
of Q 2 as soon as the photon size |r 7 |, i.e the resolution scale, becomes comparable 
to the proton size. Our model reproduces this behavior by using the perturbative 
photon wave functions with Q 2 -dependent quark masses, rrif(Q 2 ), that interpolate 
between the current (large Q 2 ) and the constituent (small Q 2 ) quark masses as 
explained in detail in Appendix |A|. The decrease of er*°* with increasing Q 2 results 
from the decreasing dipole sizes since small dipoles do not interact as strongly as 
large dipoles. 

The x-dependence of the computed proton structure function F 2 (x, Q 2 ) is shown 
in Fig. [| for Q 2 = 0.3, 2.5, 12 and 120 GeV 2 in comparison to the data measured by 
the HI [[n]] and ZEUS [|72[ detector. Within our model, the increase of F 2 (x, Q 2 ) 
towards small Bjorken x becomes stronger with increasing Q 2 in agreement with the 
trend in the HERA data. This behavior results from the fast energy growth of the 
perturbative component that becomes more important with increasing Q 2 due to 
the smaller dipole sizes involved. 

As can be seen in Fig. [8], the data show a stronger increase with decreasing x 
than the model outside the low-Q 2 region. This results from the weak energy boost 
of the non-perturbative component that dominates F 2 (x, Q 2 ) in our model. In fact, 
even for large Q 2 the non-perturbative contribution overwhelms the perturbative 
one, which explains also the overestimation of the data for x > 10~ 3 . 

This problem is typical for the SVM model applied to the scattering of a small 
size dipole off a proton. In an earlier application by Ruter |23| , an additional cut-off 



was introduced to switch from the non-perturbative to the perturbative contribution 
as soon as one of the dipoles is smaller than r cut = 0.16 fm. This yields a better 
agreement with the data also for large Q 2 but leads to a discontinuous dipole- 



proton cross section. In the model of Donnachie and Dosch |37|, a similar SVM- 
based component is used also for dipoles smaller than R c = 0.22 fm with a strong 
energy boost instead of a perturbative component. Furthermore, their SVM-based 



38 



10 5 r 



10 4 % 



Q; io 3 



10" : 



10 1 



s [GeV] 



245 (x 128) 



210 (x 64) 



170 (x 32) 



140 (x 16) 



: 115 (x 



. 95 (x 4) 



. 75 (x 2) 



2 _ 60 



■ 20 (x 0.5) 



T 1 — I 1 — r 

Photoproduction Data 

Data at y/s = 20 GeV 
ZEUS 1 



H - 



- - -o- - - 



ZEUS 2 >--^- J 
ZEUS 3 h+H 



HI 2 I — e — I 
Theory 



10" 




Q 2 [GeV' 



Figure 7: The total cross section, a l °} p {s, Q 2 ), is shown as a function of the photon 
virtuality Q 2 for cm. energies from yfs = 20 GeV to -y/i = 245 GeV, where the model 
results (solid lines) and the experimental data at different y/s values are scaled with the 
indicated factors. The low energy data at = 20 GeV are from p7|], the data at higher 
energies from the HI |38| and ZEUS collaboration [69]. The photoproduction (Q 2 = 0) 
data from p| are displayed at Q 2 = 0.012 GeV 2 . 
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Figure 8: The x-dependence of the computed proton structure function F2{x,Q 2 ) (solid 
line) is shown for Q 2 = 0.3, 2.5, 12 and 120 GeV 2 in comparison to the data measured by 
the HI p| and ZEUS p| detector, and the low energy data at yfs = 20 GeV from ||6^ 



component is tamed for large Q 2 by an additional a s (Q 2 ) factor. 

We did not follow these lines in order to keep a continuous, Q 2 -independent 
dipole-proton cross section and, therefore, cannot improve the agreement with the 
F 2 (x, Q 2 ) data without losing quality in the description of hadronic observables. 
Since our non-perturbative component relies on lattice QCD, we are more confi- 
dent in describing non-perturbative physics and, thus, put more emphasis on the 
hadronic observables. Admittedly, our model misses details of the proton structure 
that become visible with increasing Q 2 . In comparison, most other existing models 
provide neither the profile functions nor a simultaneous description of hadronic and 
7*-induced processes. 



40 



5.3 The Slope B of Elastic Forward Scattering 



The local slope of elastic scattering B(s,t) is defined as 

~do- el 



dt 



{s,t) 



(5.10) 



and, thus, characterizes the diffractive peak of the differential elastic cross section 
da el /dt(s,t) discussed below. Here, we concentrate on the slope for elastic forward 
(t = 0) scattering also called slope parameter 

m := B(M = 0, = jS<Ph\b^JMh\) = i 2 

2 /(Pftj. |6j_|) 2 W 

which measures the rms interaction radius (b 2 ) of the scattered particles, and does 
not depend on the opacity. 



We compute the slope parameter with the profile function from the T-matrix ( [2.66 ) 



and neglect the reggeon contributions, which are relevant only at small cm. energies, 
so that the same result is obtained for oh and ab scattering. 

In Fig. D, the resulting slope parameter B(s) is shown as a function of y/s for 
pp and pp scattering (solid line) and compared with the pp (open circles) and pp 
(closed circles) data from |73H75|. As expected from the opacity independence of the 
slope parameter (|5.11|) , saturation effects as seen in the total cross sections do not 



occur. Indeed, one observes an approximate B(s) oc R 2 (s) oc ln 2 ( v /i/ v /5o) growth 
for y/s > 10 4 GeV. This behavior agrees, of course, with the transverse expansion 
of J pp (s, \b±\) for increasing y/s shown in Fig. ^. Analogous results are obtained also 
for up and Kp scattering. 

For the good agreement of our model with the data, the finite width of the lon- 
gitudinal quark momentum distribution in the hadrons, i.e. Az p , Az-^, and Azk / 
in ( |A.1| ), is important as the numerator in ( |5.11| ) depends on this width. In fact, 
B{s) comes out more than 10% smaller with Az p , Az n , and Azk = 0. Furthermore, 
a strong growth of the perturbative component, e p = 0.73, is important to achieve 
the increase of B(s) for y/s > 500 GeV indicated by the data. 

It must be emphasized that only the simultaneous fit of the total cross section and 
the slope parameter provides the correct shape of the profile function J(s, \b±\). This 
shape leads then automatically to a good description of the differential elastic cross 
section da el /dt(s,t) as demonstrated below. Astonishingly, only few phenomeno- 
logical models provide a satisfactory description of both quantities |35|,|62|. In the 



approach of [14], for example, the total cross section is described correctly while the 
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Figure 9: The elastic slope parameter B(s) is shown as a function of the cm. energy 
^/s for pp and pp forward (t = 0) scattering. The solid line represents the model 
result that is compared with the data for pp (open circles) and pp (closed circles) 
reactions from [|73|-|75]|. 



slope parameter exceeds the data by more than 20% already at y/s = 23.5 GeV and, 
thus, indicates deficiencies in the form of J(s, \b±\). 



5.4 The Differential Elastic Cross Section 



The differential elastic cross section obtained from the squared absolute value of the 
T-matrix element 

da el 1 

' (s,t) = -—\T(s,t)\" (5.12) 



dt 



16tts z 



can be expressed for our purely imaginary T-matrix ( 2.66[) in terms of the profile 
function 

! [ fd 2 b ± e^ J(s,\bx\)^ 



dt M) = 4^ 



(5.13) 



and is, thus, very sensitive to the transverse extension and opacity of the scattered 
particles. Equation ( p. 13 ) reminds of optical diffraction, where J(s, \b±\) describes 
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the distribution of an absorber that causes the diffraction pattern observed for inci- 
dent plane waves. 



In Fig. ITO, the differential elastic cross section computed for pp and pp scattering 
(solid line) is shown as a function of \t\ = at yjs = 23.5, 30.7, 44.7, 63, 546, and 
1800 GeV and compared with experimental data (open circles), where the pp data 



at -y/i = 23.5, 30.7, 44.7, and 63 GeV were measured at the CERN ISR fl64f , the pp 



data at y/s = 546 GeV at the CERN SppS 0, and the pp data at -y/s = 1.8 TeV at 



the Fermilab Tevatron [75, 76]. The prediction of our model for the pp differential 



elastic cross section at the CERN LHC, \/s = 14 TeV, is given in Fig. IT 



For all energies, the model reproduces the experimentally observed diffraction 
pattern, i.e the characteristic diffraction peak at small \t\ and the dip structure 
at medium \t\. As the energy increases, also the shrinking of the diffraction peak is 
described which reflects the rise of the slope parameter B(s,t = 0) already discussed 
above. The shrinking of the diffraction peak comes along with a dip structure that 
moves towards smaller values of \t\ as the energy increases. This motion of the dip 
is also described approximately. 

The dip in the theoretical curves reflects a change of sign in the T-matrix ele- 
ment (|2.66|) . As the latter is purely imaginary, it is not surprising that there are 
deviations from the data in the dip region. Here, the real part is expected to be im- 



portant [|76| which is in the small \t\ region negligible in comparison to the imaginary 
part. 

The difference between the pp and pp data, a deep dip for pp but only a bump 
or shoulder for pp reactions, requires a C = — 1 contribution. Besides the reggeon 
contribution at small energies^ one expects here an additional perturbative C — — 1 
contribution such as three-gluon exchange |77|,[78]] or an odderon |79|-^TJ. In fact, 
allowing a finite size diquark in the (anti-)proton an odderon appears that supports 



the dip in pp but leads to the shoulder in pp reactions [31 



For the differential elastic cross section at the LHC energy, yfs = 14 TeV, the 
above findings suggest an accurate prediction in the small- \t\ region but a dip at 
a position smaller than the predicted value at \t\ ~ 0.35 GeV 2 . Our confidence in 
the validity of the model at small \t\ is supported additionally by the total cross 
section that fixes da el /dt(s,t = 0) and is in agreement with the cosmic ray data 
shown in Fig. || Concerning our prediction for the dip position, it is close to the 



value \t\ fa 0.41 GeV of [62] but significantly below the value \t\ fa 0.55 GeV of [14 . 



16 Zooming in on the result for y/s = 23.5 GeV, one finds further an underestimation of the data 
for all |t| before the dip, which is correct as it leaves room for the reggeon contribution being 
non-negligible at small energies. 
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Figure 10: The differential elastic cross section for pp and pp scattering is shown as 
a function of \t\ up to 2.5 GeV 2 . The result of our model, indicated by the solid line, 
is compared for ^/i = 23.5, 30.7, 44.7, and 63 GeV with the CERN ISR pp data @, for 
^/i = 546 GeV with the CERN SppS pp data Q, and for ^/s = 1.8 TeV with the Fermilab 
Tevatron pp data [75, 7(|, all indicated by the open circles with error bars. 
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Figure 11: The prediction of our model for the pp differential elastic cross section at 
LHC (\/s = 14 TeV) is shown as a function of momentum transfer \t\ up to 1 GeV 2 . 



Beyond the dip position, the height of the computed shoulder is always above the 
data and, thus, very likely to exceed also the LHC data. In comparison with other 
works, the height of our shoulder is similar to the one of p2j but almost one order 



of magnitude above the one of [14 



Considering Figs, pi] and |TT] more quantitatively in the small- |t| region, one can 
use the well known parametrization of the differential elastic cross section in terms 
of the slope parameter B(s) and the curvature C(s) 

da el /dt(s, t) = da el /dt(s, t = 0) exp [B(s)t + C(s)t 2 ] . (5.14) 

Using B(s) from the preceding section and assuming for the moment C(s) = 0, 
one achieves a good description at small momentum transfers and energies, which 
is consistent with the approximate Gaussian shape of J pp (s, \b±\) at small energies 
shown in Fig. 0. The dip, of course, is generated by the deviation from the Gaussian 
shape at small impact parameters. According to Q5.14 ), the shrinking of the diffrac- 



tion peak with increasing energy reflects simply the increasing interaction radius 
described by B(s). 

For small energies ^/i, our model reproduces the experimentally observed change 
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in the slope at \t\ ~ 0.25 GeV 2 [BIJ that is characterized by a positive curvature. 



For LHC, we find clearly a negative value for the curvature in agreement with | 62fl 
but in contrast to [|14]] . The change of sign in the curvature reflects the transition of 
J( s , \b±\) from the approximate Gaussian shape at low energies to the approximate 
step-function shape (|5.4|) at high energies. 



Important for the good agreement with the data is the longitudinal quark mo- 
mentum distribution in the proton. Besides the slope parameter, which characterizes 
the diffraction peak, also the dip position is very sensitive to the distribution width 
Az p , i.e. with Az p = the dip position appears at more than 10% lower values of \t\. 



In the earlier SVM approach the reproduction of the correct dip position was 
possible without the ^-dependence of the hadronic wave functions but a deviation 
from the data in the low-|t| region had to be accepted. In this low-|t| region, we 



achieved a definite improvement with the new correlation functions (|2.34j ) and the 
minimal surfaces used in our model. 

The differential elastic cross section computed for 7r ± p and K p reactions has the 
same behavior as the one for pp (pp) reactions. The only difference comes from the 
different z-distribution widths, Az n and Azk, and the smaller extension parameters, 
S n and Sk-i which shift the dip position to higher values of \t\. This is illustrated in 
Fig. |l^, where the model results (solid line) for the it^p and K^p differential elastic 
cross section as a function of \t\ are shown at yfs = 19.5 GeV in comparison with 
experimental data (closed circles) from [ 83| . The deviations from the data towards 
large \t\ leave room for odderon and reggeon contributions. Indeed, with a finite 
diquark size in the proton, an odderon occurs that improves the description of the 
data at large values of \t\ p4 



5.5 The Elastic Cross Section cr e \ a el /cr tot , and cr tot / B 

The elastic cross section obtained by integrating the differential elastic cross section 

dt^(s,t)=J (5-15) 



reduces for our purely imaginary T-matrix ( [2.66] ) to 



a e \s) = Jd 2 b ± \J(s,\b ± \)\ 2 . (5.16) 

— # 

Due to the squaring, it exhibits the saturation of J(s, \b±\) at the black disc limit 
more clearly than a tot (s). Even more transparent is the saturation in the following 
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Figure 12: The differential elastic cross section da el / dt(s , t) is shown versus the momen- 
tum transfer \t\ for ir^p and K^p reactions at the cm. energy -^/i = 19.5 GeV. The model 



results (solid line) are compared with the data (closed circles with error bars) from 82] 



ratios given here for a purely imaginary T-matrix 



" .) = S<PhMM , (5 . 17) 



° tot 2fd 2 b ± J(s,\b ± \) 
a tot (2fd 2 b ± J(s,\b ± \) 



B Jd*b ± \b ± \*J(s,\b ± \) 



(5.18) 



which are directly sensitive to the opacity of the particles. This sensitivity can be 
illustrated within the approximation 

T(s,t) = isa tot (s) exp[B(s)t/2] (5.19) 

that leads to the differential cross section (|5.14j) with C(s) = 0, i.e. an exponential 
decrease over \t\ with a slope B(s). As the purely imaginary T-matrix element ( 5.19|) 
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is equivalent to 



J(s, \b±\) = (a tot /AnB) exp[-|6 ± | 2 /2B] = (Aa el /a tot ) exp[-|6 ± | 2 /2B] , (5.20) 

one finds that the above ratios are a direct measure for the opacity at zero impact 
parameter 

J(s, \b ± \ = 0) = (a tot /AnB) = (Aa el /a tot ) . (5.21) 

For a general purely imaginary T-matrix, T(s,t) = isa tat g(\t\) with an arbitrary 
real-valued function g(|£|), J(s, \b±\ = 0) is given by (a el /a tot ) times a pure number 
which depends on the shape of #(]£])• 

We compute the elastic cross section a el and the ratios a el /a tot and a tot /B in 
our model without taking into account reggeons. In Fig. the results for pp and 
pp reactions (solid lines) are compared with the experimental data (open and closed 
circles). The data for the elastic cross section are taken from and the data for 
a tot and B from the references given in previous sections. For pp (pp) scattering, we 
indicate explicitly the prediction for LHC at yfs = 14 TeV and the onset of the black 
disc limit at yfs = 10 6 GeV. The model results for irp and Kp reactions are presented 
as dashed and dotted line, respectively. For the ratios, the asymptotic limits are 
indicated: Since the maximum opacity or black disc limit governs the yfs — ► oo 
behavior, a el /a tot (a tot /B) cannot exceed 0.5 (8ir) in hadron-hadron scattering. 

In the investigation of pp (pp) scattering, our theoretical curves confront suc- 
cessfully the experimental data for the elastic cross section and both ratios. At low 
energies, the data are underestimated as reggeon contributions are not taken into 



account. Again, the agreement is comparable to the one achieved in [62] and better 



than in the approach of 14], where a el comes out too small due to an underestima- 



tion of do Jdt in the low-|t| region. 

Concerning the energy dependence, a el shows a similar behavior as a tot but with a 
more pronounced flattening around yfs > 10 6 GeV. This flattening is even stronger 
for the ratios, drawn on a linear scale, and reflects very clearly the onset of the 
black disc limit. As expected from the simple approximation (|5.21|) , o el jo tot and 
a tot / B show a similar functional dependence on yfs. At the highest energy shown, 
yfs = 10 s GeV, both ratios are still below the indicated asymptotic limits, which 
reflects that J(s, \b±\) still deviates from the step-function shape Q5.4J) . The ratios 
a /a tot and <j tot / B reach their upper limits 0.5 and 8tt, respectively, at asymptotic 
energies, yfs — > oo, where the hadrons become infinitely large, completely black 
discs. 

Comparing the pp (pp) results with the ones for 7rp and Kp, one finds that the 



results for a tot /B converge at high energies as shown in Fig. 13. This follows from 
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Figure 13: The elastic cross section a and the ratios a /a tot and a tot /B are shown as 
a function of the cm. energy yfs. The model results for pp (pp), vrp, and Kp scattering 
are represented by the solid, dashed and dotted lines, respectively. The experimental data 
for the pp and pp reactions are indicated by the open and closed circles, respectively. The 
data for the elastic cross section are taken from [|l|] and the data for a tot and B from the 
references given in previous sections. 
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the identical normalizations of the hadron wave functions that lead to an identical 
black disc limit for hadron-hadron reactions. 



6 Conclusion 



We have developed a loop-loop correlation model combining perturbative and non- 
perturbative QCD to compute high-energy reactions of hadrons and photons. We 
have aimed at a unified description of hadron-hadron, photon-hadron, and photon- 
photon reactions involving real and virtual photons as well. Being particularly 
interested in saturation effects that manifest the S'-matrix unitarity, we have inves- 
tigated the scattering amplitudes in impact parameter space since the black disc 
limit of the profile function is the most explicit signature of unitarity. Using a lead- 
ing twist, next-to-leading order DGLAP relation, we have also estimated the impact 
parameter dependent gluon distribution of the proton xG(x, Q 2 , \b±\) to study gluon 
saturation as a manifestation of the S'-matrix unitarity at small Bjorken x. In addi- 
tion, the calculated profile functions provide an intuitive geometrical picture for the 
energy dependence of the cross sections and allow us to localize saturation effects in 
the experimental observables. 

Following the functional integral approach to high-energy scattering in the eikonal 
approximation P"T)|-|r3|l , the scattering hadrons and photons are described by light- 
like Wegner- Wilson loops (color-dipoles) with size and orientation weighted with ap- 
propriate light-cone wave functions [O. The resulting S'-matrix element factorizes 
into the universal correlation of two light-like Wegner- Wilson loops Sdd (loop-loop 
correlation function) and reaction-specific light-cone wave functions. This factor- 
ization has allowed us to provide a unified description of hadron-hadron, photon- 
hadron, and photon-photon scattering. We have used for hadrons the phenomenolog- 
ical Gaussian wave function p3, and for photons the perturbatively derived wave 



function with running quark masses rrif(Q 2 ) to account for the non-perturbative re- 
gion of low photon virtuality Q 2 



The loop- loop correlation function Sdd has been computed in the approach of 
Berger and Nachtmann ||14|| . The loop- loop correlation function Sdd has been ex- 
pressed in terms of surface (Sip) integrals over the gauge-invariant bilocal gluon field 
strength correlator. We have divided this correlator into a non-perturbative and a 
perturbative component. The stochastic vacuum model (SVM) |L5| has been used 
for the non-perturbative low frequency background field (long-distance correlations) 
and perturbative gluon exchange for the additional high frequency contributions 
(short- distance correlations) since this combination is supported by lattice investi- 



gations of the gluon field strength correlator The exponential correlation 
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function used in the non-perturbative component has been adopted directly from a 
lattice investigation of the correlator . Since this correlation function stays posi- 
tive for all Euclidean distances, it is compatible with a spectral representation of the 



correlation function [20], which means a conceptual improvement since the correla- 
tion function that has been used in earlier applications of the SVM becomes negative 
at large distances fl2|. ^4|, [2T|- f2~5[| . We have presented for the first time an explicit 



computation of the surface integrals using minimal surfaces (£1,2) bounded by the 
Wegner- Wilson loops. This surface choice is usually used to obtain Wilson's area 
law in Euclidean space |15|, 16]. Moreover, the simplicity of the minimal surfaces has 



allowed us to identify very clearly a string-string interaction in the non-perturbative 
component. 

The strongest assumption in the presented work is the form of the energy depen- 
dence introduced phenomenologically into the loop- loop correlation function Sdd- 
Motivated by the two-pomeron model of Donnachie and Landshoff 0, we have as- 
cribed to the non-perturbative and to the perturbative component a weak and a 
strong energy dependence, respectively. The constructed T-matrix element shows 
Regge behavior at moderately high energies and contains multiple gluonic interac- 
tions important to respect unitarity in impact parameter space at ultra-high ener- 
gies. 

The model parameters have been adjusted to reproduce a wealth of experimen- 
tal data (including cosmic ray data) over a large range of cm. energies: total, 
differential, and elastic cross sections, structure functions, and slope parameters - 
including cosmic ray data. The model parameters that allowed a good fit to high- 
energy scattering data are in good agreement with complementary investigations: 
The parameters of the non-perturbative component — the correlation length a, the 
non-Abelian strength k, and the gluon condensate G<i — are constrained by lat- 
tice QCD investigations, by the string tension o of a static quark-antiquark pair, 
and by the SVZ gluon condensate G2 essential in QCD sum rule investigations. 
The parameters of the perturbative component have not been adjusted. We have 
used the p-meson mass for the effective gluon mass representing the IR regula- 
tor and have determined M 2 so that the strong coupling freezes for large distance 
scales at a s = 0.4, where our non-perturbative component is at work according to 
a low energy theorem [16|,|33|]. For the energy dependence, the exponents of the 



Donnachie-Landshoff two-pomeron fit, e so ft and ehard, have been used as an orien- 
tation for our energy exponents e NP and e p . Besides these parameters describing 
the universal loop-loop correlation function Sdd-, the reaction-dependent parameters 
in the light-cone wave functions are also consistent with other approaches: In the 
hadron wave functions, the transverse extension parameters Sh have been found in 



good agreement with the corresponding electromagnetic radii [25[ while the width of 
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the longitudinal quark momentum distributions Azh has been computed from ]zb 
In the photon wave function, the running quark masses that coincide with the cur- 
rent quark masses for large Q 2 and the constituent quark masses for small Q 2 |22j 
have been chosen in agreement with sum rule derivations. 

Having adjusted the model parameters, we have studied S'-matrix unitarity limits 
of the scattering amplitudes in impact parameter space. On the basis of dipole-dipole 
scattering, we have found explicitly that our model preserves the unitarity condition 
and attains the black disc limit at ultra-high energies. The profile functions have 
been calculated for proton-proton and photon-proton scattering and have shown 
very clearly that the interacting particles become blacker and larger with increasing 
energy. At ultra-high energies, the opacity saturates at the black disc limit while the 
transverse expansion of the scattered particles continues. Moreover, in longitudinal 
photon-proton scattering, we have observed that with increasing photon virtuality 
Q 2 not only the maximum opacity increases but also the energy at which it is reached 
for zero impact parameter. 

Using a leading twist, next-to-leading order QCD relation between the gluon 
distribution of the proton xG(x, Q 2 ) and the longitudinal proton structure function 
Fl(x,Q 2 ) P8|,fl9|l, we have related the impact parameter dependent gluon distribu- 
tion xG(x, Q 2 , \b±\) to the profile function for longitudinal photon-proton scattering 
and found low-x saturation of xG(x,Q 2 , \b±\) as a manifestation of S'-matrix uni- 
tarity. In accordance with the profile function, xG(x,Q 2 , \b±\) decreases from the 
center towards the periphery of the proton. With increasing photon virtuality Q 2 , 
the increase of xG(x, Q 2 , \b±\ = 0) becomes stronger towards small x and the satu- 
ration value of xG(x,Q 2 , \b±\ = 0) increases but is reached at decreasing values of 
x. In contrast, at fixed Q 2 , the integrated gluon distribution xG(x,Q 2 ) does not 
saturate because of the growth of the proton radius with decreasing x observed in 
our approach. Similar results are obtained in complementary approaches H^7|, |50|, |51 



More model dependent are the specific energies where these saturation effects 
set in. The profile function saturates the black disc limit at zero impact parameter 
for y's > 10 6 GeV in proton-proton scattering and for > 10 7 GeV in longitudinal 
photon-proton scattering with Q 2 > 1 GeV 2 . In both reactions, the wave function 
normalization determines the maximum opacity. The saturation of xG(x,Q 2 , \b±\) 
occurs in our approach for Q 2 > 1 GeV 2 at values of x < 10 -10 , far below the HERA 
and THERA range. 

For proton-proton scattering, we have observed that the rise of the total and 
elastic cross section becomes weaker for ^Js > 10 6 GeV due to the onset of the 
black disc limit at = in the profile function. This saturation of the profile 
function has become even more apparent in the ratios o el /a tot and a tot / B which are 
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a measure of the proton opacity. In contrast, no saturation effect has been observed 
in the slope parameter B(s) which is a measure of the transverse expansion of the 
proton. Considering the differential elastic cross section da el /dt, the model has 
described the diffraction pattern and also the shrinkage of the diffraction peak with 
increasing energy in good agreement with experimental data at small momentum 
transfers \t\. Around the dip region, where a real part is expected to be important, 
deviations from the data have reflected that our T-matrix is purely imaginary. Our 
predictions for proton-proton scattering at LHC (y/s = 14TeV) have been a total 
cross section of cr*°* = 114.2 mb in good agreement with the cosmic ray data and a 
differential elastic cross section da el /dt with a slope parameter of B = 21.26 GeV~ 2 , 
a negative curvature, C < 0, and a dip at \t\ ~ 0.35 GeV 2 . 

For pion-proton and kaon-proton scattering, results analogous to proton-proton 
scattering have been obtained but with a slightly stronger rise observed in the total 
cross section. This has been traced back to the smaller size of pions and kaons in 
comparison to protons, S p = 0.86 fm > S w = 0.607 fm > Sk — 0.55 fm, and the 
perturbative component becoming increasingly important with decreasing dipole 
sizes involved. Furthermore, a weak convergence of the ratios cr^/a^ and 0"*°* /cr^ 
towards unity has been observed as the energy increases. The smaller size of the 
pion and kaon has also been reflected in the differential elastic cross sections da el /dt, 
where the dip is shifted towards larger values of \t\. 

For photon-proton and photon-photon reactions, an even stronger rise of the 
total cross section has been observed with increasing energy. As illustrated in the 
proton structure function F 2 (x, Q 2 ), this rise becomes steeper with increasing photon 
virtuality Q 2 . Again, we have traced back the strong energy boost to the growing 
importance of the perturbative component with decreasing dipole size. Besides some 
deviations from the experimental data with increasing Q 2 , our model has described 
cr*°p(s, Q 2 ) successfully in the low-Q 2 region where the running quark masses become 
constituent quark masses. 

We plan to present complementary investigations within our model in future 
work. Going to momentum space, we calculate on the two-gluon exchange level 
the unintegrated gluon distribution of the proton. Insight into the non-perturbative 
structure of this distribution can be gained from the non-perturbative component 
of our model [[H]]. In Euclidean space-time, we compute the static quark- ant iquark 



potential, the associated flux tube, and the van der Waals force between two static 
color-dipoles flB]] . It is a long-range project to implement the energy dependence 



more fundamentally. In a recent attempt, the energy dependence of high-energy 
scattering has been related successfully to critical properties of an effective near 



light-cone Hamiltonian in a non-perturbative lattice approach [85]. Furthermore, the 



correlation of inclined Euclidean Wegner- Wilson loops generates energy dependence 
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after analytic continuation to Minkowski space-time ||86|| . Here, encouraging new 
results obtained with instantons and within the AdS/CFT correspondence | 88[ 
have to be compared to calculations in the stochastic vacuum model. 
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A Hadron and Photon Wave Functions 



The light-cone wave functions ^(z^rj) provide the distribution of transverse size 
and orientation r*j and longitudinal quark momentum fraction Zi to the light-like 
Wegner- Wilson loops W[Cj] that represent the scattering color-dipoles. In this way, 
they specify the projectiles as mesons, baryons described as quark-diquark systems, 
or photons that fluctuate into a quark-antiquark pair before the interaction. 



The Hadron Wave Function 



In this work, mesons and baryons are assumed to have a quark-antiquark and 
quark-diquark valence structure, respectively. As quark-diquark systems are equiv- 
alent to quark-antiquark systems this allows us to model not only mesons but 
also baryons as color-dipoles represented by Wegner- Wilson loops. To characterize 
mesons and baryons, we use the phenomenological Gaussian Wirbel-Stech-Bauer 



ansatz [26 | 



iph(zi,ri) 



Zj(l- Zi) 

2nS 2 h N h 



D 2 /(4A4) , 



\fi\ 2 /(iS 2 h ) 



(A.l) 



where the hadron wave function normalization to unity 



dzid 2 ri \ipi(zi,fi)\ : 



(A.2) 
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Table 1: Hadron Parameters 



Hadron 


Az h 


S h [fm] 


P,P 


0.3 


0.86 




2 


0.607 


K± 


0.57 


0.55 



requires the normalization constant 

N h = C d Zi Zi(l - e -(^-D 2 /(2A^) _ (A _ 3) 



The different hadrons considered — protons, pions, and kaons — are specified by 
Azh and Sh providing the width for the distributions of the longitudinal momentum 
fraction carried by the quark Z( and transverse spatial extension \fi\, respectively. 
In this work, the extension parameter Sh is a fit parameter that should resemble 



approximately the electromagnetic radius of the corresponding hadron | 2~5| , while 
Azh = io/(v2w/j) |26|] is fixed by the hadron mass vrth and the value w = 0.35 — 
0.5 GeV extracted from experimental data. We find for (anti-)protons Az p = 0.3 
and S p = 0.86 fm, for pions Az n = 2 and S n = 0.607 fm, and for kaons Azk = 0.57 
and Sk = 0.55 fm which are the values used in the main text. For convenience they 
are summarized in Table 0. 

Concerning the quark-diquark structure of the baryons, the more conventional 
three-quark structure of a baryon would complicate the model significantly but 
would lead to similar predictions once the model parameters are readjusted (12 



In fact, there are also physical arguments that favor the quark-diquark structure of 
the baryon such as the 51 = 1/2 enhancement in semi-leptonic decays of baryons f8*9 



and the strong attraction in the scalar diquark channel in the instanton vacuum ]90| . 



The Photon Wave Function 



The photon wave function ^(zi, r*j, Q 2 ) describes the fluctuation of a photon with 
virtuality Q 2 into a quark-antiquark pair with longitudinal quark momentum frac- 
tion Zi and spatial transverse size and orientation r*j. The computation of the corre- 
sponding transition amplitude (qq{zi,fi)\ , y*(Q 2 )) can be performed conveniently in 
light-cone perturbation theory |HJ and leads to the following squared wave functions 
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for transverse (T) and longitudinally (L) polarized photons |34| 



\^ T { Zi A,Q 2 )? = { K + - ^) 2 ] 4^/ 1*1) + ™/#o(*/ lnl)}(A-4) 

hM*. r i5 Q 2 )| 2 = |^E e ?( 4 ^ 2 * " 1*1)} , (A.5) 

where a is the fine-structure constant, ej is the electric charge of the quark with 
flavor /, and Kq and K\ are the modified Bessel functions (McDonald functions). 
In the above expressions, 

e) = z t (\ - Zi ) Q 2 + m) (A.6) 

controlls the transverse size (-distribution) of the emerging dipole, \fi\ oc 1/e/, that 
depends on the quark flavor through the current quark mass rrif. 



For small Q 2 , the perturbatively derived wave functions, ( |A.4| ) and ( |A.5[ ), are 



not appropriate since the resulting large color-dipoles, i.e. \fi\ oc 1/rrif ^> 1 fm, 
should encounter non-perturbative effects such as confinement and chiral symmetry 
breaking. To take these effects into account the vector meson dominance (VMD) 
model |^2j is usually used. However, the transition from the "partonic" behavior at 
large Q 2 to the "hadronic" one at small Q 2 can be modelled as well by introducing 
Q 2 -dependent quark masses, m/ = rrif(Q 2 ), that interpolate between the current 
quarks at large Q 2 and the constituent quarks at small Q 2 ||22|| . Following this 
approach, we use (|A.4j) and (|A.5|) also in the low-Q 2 region but with the running 



quark masses 

muAQ 2 ) = 0.178 GeV (1 - 0(0^ - Q 2 ) , (A.7) 

m s (Q 2 ) = 0.121 GeV + 0.129 GeV (1 - ^) Q(Q 2 S - Q 2 ) , (A.8) 
and the fixed charm quark mass 

m c = 1.25 GeV , (A.9) 

where the parameters Q 2 ud = 1.05 GeV 2 and Q 2 = 1.6 GeV 2 are taken directly 
from [^2[ while we reduced the values for the constituent quark masses rrif(Q 2 = 0) 
of |22] by about 20%. The smaller constituent quark masses are necessary in order 



to reproduce the total cross sections for 7*p and 7*7* reactions at low Q 2 . Similar 
running quark masses are obtained in a QCD-motivated model of the spontaneous 



chiral symmetry breaking in the instanton vacuum |93[ that improve the description 
of 7*p scattering at low Q 2 |94 . 
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B Correlation Functions 



In this Appendix, we describe explicitly the way from the simple exponential correla- 
tion functions in Euclidean space-time ( |2.34| ) to their transverse Fourier transforms 
in Minkowski space-time, ( [2.50|) and ( [2.55|) . The first step in this procedure is 
the Fourier transformation of the exponential correlation functions ( |2.34| ) in four- 
dimensional Euclidean space-time 

D E {K 2 ) = D E {K 2 ) = Jd 4 ZD E {Z 2 /a 2 )e lKZ 

POO PIT PIT /*27T 

= / d\Z\ \Z\ 3 / d(j) 3 sin 2 3 / d(j) 2 sin0 2 / dfa D E {Z 2 /a 2 ) e -W*l«**s 
Jo Jo Jo Jo 

r d\Z\\Z\ 2 D E (Z 2 /a 2 )J 1 (\K\\Z\) = - l 271 " . , (B.l) 



K\ Jo a [K 2 + a~ 



2 



where J% is the 1st order Bessel function of the first kind. Here, the Euclidean metric 
—8nv and four-dimensional polar coordinates and the corresponding four-volume 
element d A Z = d\Z\ \Z\ 3 d(f>3 sin 2 03<i02 sin02 d(j>\ have been used. With (P-l| ), one 
obtains 

= (B „ 

Now, (p.l| ) and ( p.2|) are analytically continued to Minkowski space-time, K± — > ik° 
or equivalently -8^ -> = diag(l, -1,-1,-1), 

D(k*) = 1 ^— T , DUP) = . (B.3) 

a (— k 2 + a 2 )z a(— k 2 + a 2 )2 

Setting k° = k 3 = 0, which is enforced in the computation of \ by 5-functions, one 
finds k 2 = —k\ and consequently 

D^Ck 2 ± ) = , bT\kl) = - 7 . (B.4) 

a (k 2 + a -2 ) 2 a(/c 2 + a -2 ) 2 



The transverse Fourier transformation (|2.48|) of these two expressions is the remain- 



ing step that leads directly to (|2.50|) and (|2.55|). 
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